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SYNOPSIS 

qnWlE ASPECTS OP QUANTUM THEORY OE CHANNELING 
Anand Prakash Pathak 
Ph.D. 

Department of Physics 
Indian Institute of Technology Kanpur 
March 1971 

In the present work, the quantum mechanical approach to the 
theoiy of channeling and blocking of particles in crystals has 
been investigated using Debye model for lattice vibrations and 
one phonon approximation for inelastic processes taking place in 
the crystal when the particles propagate through it. The channel- 
ing and emission of charged particles in crystals and anomalous 
penetration of neutrons have been considered to illustrate the 
various aspects of the quantum theoiy of channeling. The rate of 
energy loss of charged particles due to conduction electrons is 
important under channeling conditions and has been calculated 
using a dielectric formulation. 

After a review of the classical theoiy of Lindhard in the 
Chapter I, the need for a quantum theoiy has been briefly dis- 
cussed. m the Chapter II., a general formulation for channeling 
of particles into crystals has been developed and an expression 
for the renormalization of the initial state wave function due 
to inelastic processes has been obtained using a general inter- 
action potential. The approximation involved here is that the 
crystal is initially in a veiy low lying state and only the states 
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close to the initial state of the ciystal are excited hy the inci- 
dent particle trader channeling conditions. The general formtala 
has been simplified using first Bom approximation and Behye model 
for lattice vibrations to calculate the imaginary part of the renon 
malization matrix which governs the attenuation of the particles. 


In Chapter III, different potential models have been employe' 
to illustrate the details of the calculation. Using a potential 
for neutrons which goes over to Beimi pseudo-potential as its width 
vanishes, it has been shown that for anomalous penetration, the 
interaction potential should be weak compared to the particle energj 
and should be localized around the lattice sites. It is shown that 
the more localized and weak is the interaction potential, the more | 
pronounced is the channeling effect. Bor charged particles, two | 
potential models have been employed, the screened Coulomb and Bom-^ 
Mayer. The channeling behaviour is again determined by the range 
of the potential. "When the width of screened Coulomb potential 
is equal to the width of Bom-Mayer potential, the screened Coulomb 
potential is more favourable for channeling. The tempemture de- 
pendence of the process due to lattice vibrations, comes primarily 
through the Debye-Waller factor and tends to reduce the anomalous 
efffect. The energy dependence of the expressions supports the 
e XI 13 rime nt ally established fact that thickness into the 

ciMafal at which one half of the initially channeled particles 
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The problem of emission of particles from the source embedded 
in the crystal has been studied quant-um mechanically in Chapter IV. 
The real part of the renormalization matrix adds a correction term 
to the Fourier transform of the potential, which determines the 
width of the emission pattern. The magnitude of this correction 
teim is extremely small at high particle energies but starts increa. 
ing as 1/E for low energies. Unlike the previous quantum mechani- 

ir 

cal treatment of this problem, the attempt here has been to calcu- 
late the attenuation due to inelastic processes using renormaliza- 
tion matrix and is found to be small compared to phenomenologically 
chosen value. The attenuation depends upon the particle mass m^ 
and energy as (m^/E ^) An interesting feature of the pre- 
sent formalism is qualitative indication of difference in widths 
of electron and positron emission patterns even in the two beam ; 
theory. | 

I 

The energy loss of the particles moving in channels is \ 

mainly through the excitation of the conduction electrons as collecj 
tive oscillations. The corresponding contribution has been calcu- 
lated numerically in Chapter IV using a dielectric formulation 
with different models for dielectric function. The stopping power 
dE/dx is found to be very nearly inversely proportional to velocity 
for particle velocities less than about ten times the Eermi velo- 
city. For higher velocities, the usual velocity dependence 
(dE/dx cc 1/V ) is obtained. Finally the conclusions have been 
summarized in the last chapter. 


The problem of emission of particles from the soxirce emhedde 
in the ciystal has been studied quantum mechanically in Chapter IV. 
The real part of the renormalization matrix adds a correction term 
to the I'ourier transform of the potential, which determines the 
width of the emission pattern. The magnitude of this correction 
term is extremely small at high particle energies but starts increa 
ing as 1/E for low energies. Unlike the previous quantum mechani- 

ir 

cal treatment of this problem, the attempt here has been to calcu- 
late the attenuation due to inelastic processes using renormaliza- 
tion matrix and is found to be small compared to phenomenologically' 
chosen value. The attenuation depends upon the particle mass m^ ; 
and energy Ep as (m^/E^)**/^, In interesting feature of the pre- ' 
sent formalism is qualitative indication of difference in widths 
of electron and positron emission patterns even in the two beam i 
theory. 

I 

The energy loss of the particles moving in channels is 
mainly through the excitation of the conduction electrons as collec- 
tive oscillations. The corresponding contribution has been calcu- [ 

i 

lated numerically in Chapter IV using a dielectric formulation I 

I 

with different models for dielectric function. The stopping power | 
dE/dx is found to be very nearly inversely proportional to velocity! 
for particle velocities less than about ten times the Eermi velo- 
city, For higher velocities, the usual velocity dependence 
(dE/dx oc 1/V ) is obtained. Finally the conclusions have been 
summarized in the last chapter. 
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INTRODUCTION 


1.1 Discove:^ and Applications of Channeling: 

The study of channeling and blocking phenomena in ciystals 
is a new and developing branch of Solid State Physics. When 
energetic particles incident on a ciystal move along directions 
close to the principal crystallographic axes or planes, they 
penetrate anomalously long distances and are said to be channeled. 
The importance and applications of channeling technique in Solid 
State Physics is attracting many workers to look into the details 
of the phenomena and to search for new applications of the tech- 
nique. A substantial amount of experimental work has been done 
during the last five years; mostly at the Chalk River Nuclear 
Laboratories, Canada and the University of Aarhus, Denmark. o?he 
technique is finding extensive use and has proved to be very 
useful in Nuclear Physics too, as for example in nuclear life- 
time measurements. 

Historically, the directional effects were anticipated 

■1 

as early as 1912 by Stark and Wendt but the detailed theoretical 
and experimental investigation of the problem was delayed until 

p 

the last decade. In 1963, Robinson and Oen made a ccmiputer 
simxilation of the heavy ion motion in crystals. They investigated 
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the slowing down ol 1-10 keY heawy ions in a lattice model which 
included realistic repulsive interaction and obtained veiy large 
penetrations for particles with initial velocities nearly parallel 
to close-packed directions. Experimentally this effect was first 
observed in studies on the penetration of keY ion in crystals by 
Piercy et al. ^ and by Lutz and Sizemann^. They found by measur- 
ing the distribution of penetration distances, that a significant 
fraction of ions incident along low index directions had anoma- 
lously long range. At about the same time, Nelson and Thompson^ 
observed anomalously high transmission of 50-75 ke7 protons and 
helium ions in low index directions of thin gold single crystals. 
The effect was further confirmed by the experiments of Deamaley 
and others” who reported that in the transmission experiments, 
using thin single ciystals of Si, 1-10 MeY protons and alpha 
particles lose energy more slowly, if injected parallel to a 
close-packed axis or plane. Similar reductions in energy loss 

n 

weie observed' for bjcomine and iodine ion beams at energies 
upto 100 MeY. These experiments and the detailed range measure- 

Q 

ments of Komelson et al. on heavy ion beams in single crystals 
fully confirmed the channeling concept and established that in 
an aligned ciystal, more than 50% of the beam becomes channeled 
and penetrates much further than the observed range for an irre- 
gular array of target atcxQs (i.e. amorphous target). 

The first accurate and comprehensive theoretical treat- 
ment of these directional effects was given by Lindhard^ in 1965* 
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This ta:?eatment is "based on the classical orbital picture of the 
actual steering mechanism for the motion of the particle in a 
crystal. The experimental development was extremely rapid follow- 
ing this classical theory of lindhard. A good amoimt of experi- 
mental work has been done on heavy ion ranges in tungsten single 
1 0 

crystal , wide angle Rutherford scattering in tungsten and sili- 

11 12 
con t and other diamond-type lattices (such as diamond, Si, 

Ge, GaP, GaAs, GaSb), x-ray production in A1 and Cu by proton 
bombardment and neutron production by nuclear reaction of protons 
in copper ^ (Gu ^ (p, n) Zn ^). A more detailed account of expe- 
rimental work has been reported at the International Conference 

on Electromagnetic Isotope Separators and related ion Accelera- 

15a 

tors and their Application to Physics and at the International 
Conference on Atomic Collisions and Penetration Studies with 
energetic (keT) ion beams held at Chalk River, Ontaxio”*^^. 

The phenomenon of ’channeling' has proved to be a very 

useful tool for making some very important experimental studies. 

It can be used to find out whether a doped impurity has occupied 

1 17 

a substitutional site or an interstitial site ’ ' in the crystal. 

1 P> 

Further, one can study the lattice disorders and surface 
19 

effects using this technique. It has been applied to the study 

20 

of nuclear-life times also, in a very interesting way. One 
particularly useful application of the channeling technique has 
been the study of ion implantation in semiconductors^ ’ * ^ in 

which the location of the implanted ions and a determination of 
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damage produced is studied "by means of the channeling hehaviour 

of MeV projectiles. The technique has also been used in finding 

24 - 

out the location of inert gas atoms in some ciystals (such as 

KCl, Cal’2 ^‘3. UO2). Another related phenomenon, which on one 

hand, gave rise to the feeling of the existence of the directional 

effects in crystals and initiated its discovery and now, on 

the other hand, uses channeling as a technique for its own detailed 

27 

Investigation and development , is the phenomenon of ’sputtering’ 

where the yield of the atoms from the target due to homhardment 

of high ener^r ions and protons depends upon the crystallographic 

directions. A simple theory of ’sputtering’ has been given by 
28 

Lehman . A more detailed information regarding the discovery 

and applications of channeling and blocking phenomena are dis- 

29 

cussed in the review articles by Katz et al. ^ , Erginsoy, and 
30 

Davies . 

The. following section contains a brief review of the clas- 
sical theory proposed by Lindhard. The underlying basic assump- 
tions which lead to a consistent approximation method for treat- 
ing the channeling phenomena, have been given and the conditions 
of their validity are derived briefly. The applicability of 
this theory and the cases,v^here a quantum mechanical theory has 
to be used in view of the fact that the former cannot give any 
thing beyond the gross features of the phenomena, are indicated 
in the last section of the present chapter. 


0 
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1.2 Rev iew of the Classic^ fheoiy : 

The classical theory of directional effects in perifect 

q 

crystals proposed hy Lindhard^ is based on the 'orbital picture’ 
of fast moving particles and their steering by the densely 
packed strings or planes of atoms existing in the crystal. The 
theory has been constructed on the basis of four main assiamp- 
tions which lead to a consistent approximation procedure for the 
treatment of motion of the particle in channels. 

Rirst, the angles of scattering of the particle may be 
assumed to be small. This assumption is valid for fast heavy 
ions moving close to the crystallographic directions and obviously, 
scattering by large angles would imply that the original direc- 
tion is completely lost. The deflections can then be computed 
classically as if the atoms ere infinitely heavy i.e. as if the 
particles were moving in an external fixed potential. 

The second assumption is that the strongest correlation 
occurs for a collision with a row of atoms. The other atoms in 
the crystal have a negligible influence on the correlations 
occuring in one particular row of atoms. This assumption holds 
because a collision demands that the particle comes close to the 
atom and collisions are strongly correlated if the particle moves 
at a small angle with a row of atoms (as a consequence of the 
first assumption). However, if it passes close to one atom in the 
row, it must also pass close to the neighhouring atoms in the 
same row. This leads to the concept of a string of atoms. 
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cTiaracterized merely by tbe constant distance of separation 'd* 
of atoms, placed on a straight line. In first approximation, it 
can be said that the collisions occur with one string at a time, 
the string collisions being independent and uncorrelated. The 
physical importance of the string is emphasized by the fact that 
practically all the physical processes caused by the particle or 
influencing its path, demand that it comes close to the string. 

One exception is the resonance excitation of the atomic electrons 
which may take place far away from the particle, if it has a 
high velocity. The simplicity of the string approximation is 
due to the fact that the lattice structure does not enter into 
the problem, the only lattice parameter coming into picture 
being the distance d between atoms in the string. The strings 
belonging to low index directions have a small value of d, and 
are the most pronounced ones. Correlations weaker than those 
of strings are expected for crystal planes, atomic pairs etc. 

The third assumption is regarding the adequacy of the 
classical orbital picture. Since the individual collisions of 
the particles with charge z^e and velocity v with the crystal- 
atoms of charge Z2e, need quantal corrections either when the 
quantity x = (2z^Z2e /hv) is not large compared to unity or 

when the impact parameter is large? the classical approximation 
might seem to be doubtful in several cases. However, it has 
been shown by Lindhard (Bef.Sj Appendix B) that the classical 
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description of many successive collisions with atoms in a string 
does not become invalid at high velocities. 

Lastly, the idealized case of a perfect lattice and a 
perfect string may be used as a first approximation. The ther- 
mal and zero point vibrations of atoms can be considered to be 
comparatively small and their influence on the symmetry proper- 
ties of the lattice is assumed to be negligible. 

With these assumptions Lindhard introduced the idea of 
continuum approximation to the potential of a string. Using this 
in conjianction with the assumption of small angle scattering 
leads to an expression for critical angle for channeling. The 
basis of the continuum approximation is to assume that many con- 
secutive atoms in the string simultaneously influence the parti- 
cle trajectory. The average potential responsible for the deflec- 
tion is obtained by averaging over the whole string. Thus if we 
take ions as the incident particles, the average potential at a 
distance r from the string is given by 

U(r) = / ^ V (/z^ + r^') (1.1) 


where V(R) is the ion-atom potential and d is the distance between 
atoms in the string which has been assumed to be along the z-azis. 
Ror R not very much larger than a, the potential V(R) is essen- 
tially of Thomas Fermi type and we may write 

.2 


?(R) = 


z^Z2e 


4 o (R/a) 


( 1 . 2 ) 
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where a is the screening length of the particle atom interaction 

and is gi'ven by a = x 0.8853 and ^^(R/a) 

32 

is the Fermi function belonging to ..one isolated atom . Equation 
(1.1) shows that, the variation of U with 1/r is by one power 
less than the variation of V with 1/r. From the equations (1.1) 
and (1.2), we may write 



U(r) = j ^ (r/a) 

(1.3) 

with 

OO -y..' 

c (r/a)= / ^ M z^+rVa) 

(1-4) 


A detailed and more accurate estimate of C (r/a) is given 
el^where^^. For the present purpose, we may use some what simp- 
ler estimates. In order to get qualitative insight in the beha- 
viour of C, i.e. of U, we note that ^ ^ R-/^- ”1 so 

that from (1.4), ^(r/a) should increase logarithmically for 
small r. Ihus, 

C (r/a) % 2 log ~ for r < Ca (1-5) 

where the constant of integration 2 log G is determined by the 
screening. A better estimate valid for all r has been taien to be 

? (r/a) % log r + 1 3 (1.6) 

so that 5 (r/a) ^ (Ca/r) for r > Ca. The constant C has been 
chosen to be ■f3 which gives fairly good overall fit. The expres- 
sions for the density of electrons p (R) and atomic potential 
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V(R) corresponding to formula (2.6) for 5 (r/a) are given by 

,2 


P (R) = -2- s/2 

4-n; 2 ^ j5/2 


V(R) = ;2-r 




( 1 . 6a) 

(1.6b) 


The condition for validity of continuum approximation for 
the potential due to the string is obtained by demanding that 
the distance d between successive atoms in a string is small 
compaired to At. vy where Vjj = v cos is the velocity component 
paraillel to the string and can be taken approximately equal to 
V for small ^ , and collision time At = ^min^^^/^^ ) s 1 

being the impact parameter with the string and (l),the cor- 

responding minimum distance of approach. The condition for vali- 
dity of the continuum approximation thus becomes 

r_.. (1) 

At. V cos Tp ^ — — - > d (1.7) 

let us apply the condition (1.7) in its most restrictive form 
so that we demand its fulfilment for 1=0 and determine 
^min (hereafter written simply as hy 

U^^min) = 2' M^v^sin^^, f ^ (1.8) 

so that the condition (1.7) becomes 

^ exp ( - ) > 1 (1,9) 

2 

with b = z^Z 2 ® particle energy. For increasing 
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from zero, the inequality is ■violated first by the rapid decrease 
of the exponential, provided (Ca/’i'd) can remain large. Thus the 
condition (1.9) requires that 

^ = ]l2h/E^ = = 2z^Z2eVc1 (1.10) 

provided Ca/ii'-]d > 1 i.e. Tp ^ < a/d or > S' = (1.11) 

At lov/ energies whei?e (1.10) is not valid, one uses (1.6) 
in (1.7) and ( 1 .8), getting the condition 

i|; < = (Cai|;^/df2)''/^ (1.12) 

and since C/ir2 Y’3/V2 ^ 1, the critical angle ip 2 applies when 

> a/d or E < E* (1.13) 

From the above expressions (1.8)-( 1. 13) , we note that the poten- 

tial energy barrier E t-j (for E > E’) is independent of energy 

2 1/2 

but for Ep < E’, the corresponding barrier Bp’i '2 "Varies as E^ ' 
and decreases as the particle energy decreases. In contrast to 
in (1.10), ^2 (1.12) depends upon the atomic radius and 

on the behaviour of screened atomic potential. As a consequence, 
(1.12) cannot be expected to hold accurately at very low energies. 

The behaviour of a beam of particles moving through a 
lattice is qualitatively discussed by equations (1.10) and (1.12), 
If the initial angle t is less than the continuum picture 

of the string applies' and the beam obeying this condition is said 
to be aligned beam (the condition for aligned beam being t 
G’ of order 1-2) when f > C’t , one ^ts the random beam. Thus 
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the yield of a reaction requiring the incident particle to come 
close to the lattice atoms (within about 0.2 S) such as Ruther- 
ford scattering, nuclear reactions etc., will depend upon the 
angle of incidence. It starts decreasing when the angle of 
incidence with respect to a prominent ciyst alio graphic direc- 
tion is less than ^ and it reaches minim-urn -value for = 0 , 
corresponding to a completely aligned beam. This minimum yield 
has been estimated to be 


N du ( p ^ ) 

x 


where is the mean square amplitude of lattice "vibrations 


perpendicular to the string and a is the Thomas Rermi screening 
distance. 


Another important case of channeling is the case of planar 
channeling which occurs when the particle moves parallel to a 
crystallographic plane. The correlations here are weaker and 
less effective than for the string case. Using continuum appro- 
ximation, the average potential in this case can be written as 

°° '00 

Y(y) = Nd / 2Kr dr V(f ) (1.14) 

^ o 


whei?e y is the distance from the plane and Nd represents the 

Jr 

average nimber of atoms per unit area of the plane, d_ being 
the distance between the consecutive planes. The ion-atom 
potential ¥(R) is given by (1,2), Using the standard atomic 
potential (1.6b), one gets 
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1 /2 

Y(y) = 2itz^Z2e2 [(y^ + cV) - y ] (1.15) 

which corresponds "bo the expression (1«6) for ^ (r/a) for the 
string case. 


The criteria for the use of continuum potential for planar 
case is much more involved than for a string. The approximate 
condition for this case, as found hy Lindhard^, is 




-? ] > 


a 


With a = Bp /Ep, ; = y^i^/Oa, 


E ’ = 


2^226 


'P 2TtC^a^Nd, 


( 1 . 16 ) 


-V Z^Z 2 X 30 eT. 


and 


^min minimum distance of approach to the plane. 


Another consequence of directional effects is observed 
when the particles, emitted from sources (such as radioactive 
nuclei) embedded in a crystal propagate through it. The inten- 
sity of the emission pattern is found to depend strongly on the 
direction of observation with respect to a prominent crystallo- 
graphic direction or plane. For positrons, the emission pattern 
has a dip along the directions along which maximum channeling 
occiu?s (when the emitter is at a lattice site). On the other 
hand, one gets a strong peak along these directions, for the 
case of electron emission. This phenomenon is called blocking 
and has also been treated by the classical theory of lindhard. 
Supposing that the emitting nucleus is in the neighbourhood of 
an atomic position in the perfect string, one takes the probabi- 
lity distribution dp(r) of this nucleus, in space to be of the 
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G-aussian type 

2 2 

dP = e"^ 2rdr. a/p^ (1.17) 

2 

where p is the mean square displacement due to the thermal 

^2 ^ ^ 'I 

vibrations and a=(1-e ° ) is normalization constant 

and r^ is defined as mr^ = (IId)“ . Ihen the probability dis- 
tribution of transverse energy of particles is found to be 


(assuming P << r^ i.e. a % 1) 


n 


= exp (e 




-1) J-expJ^- ('^* 


and the integrated dip in the emitted intensity is given by 

2 2 2 2 
fC a. + p^ 

12 jv 2" "■ 2 ’ (1.19) 

This expression shows that as the thermal vibrations tend to 

2 

increase the mean square displacement p , the dip becomes less 
pronounced. Por the ideal case of no thermal vibrations, one 
gets maximum dip. 


Equation (1.19) can be used to calculate ^ if P is known. 
On the other hand if ii is experimentally measured then the cal- 
culated value of p from eqn. (1.19) can be compared with its 
value determined otherwise. Apart from such estimates, it is 
seen that is proportional to = 2z^Z2e /dS. Thus the par- 

ticles emitted with less energy will show a more paxjnounced dip 
than that exhibited by particles emitted with hi^er enexgy. 


18) 


14 


1,3 Need for Quotum Theo2^_? 

The classical theory reviewed in the last section as pro- 
posed hy Lindhard has been found to be in good agreement with 
experiments on heavy ion and alpha particle channeling. The mag- 
nitude of critical angle, 2?anges of ions, yields of Rutherford 
scattering and nuclear reactions etc, have been measured and 
results compared with the classical predictions a 2 ?e found to be 
in satisfactory agreement " However, for light particles, 
such as the electrons and positrons, the classical theory cannot 
be used if one needs anything beyond the gross features, mainly 
because of two reasons: firstly, the phenomena must be closely 
coupled to the Bragg reflections and other related interference 
effects which are obseinred in electron micrographs. Actually the 
phenomena observed in the electron microscope is described in 
temns of wave interference; Bragg angles and resonance widths 
dominate the intensity patterns even when the electron wave-length 
is very much smaller than the lattice spacing of the target crys- 
tal. Secondly, according to classical mechanics, the trajectory 
of a particle with a definite energy, moving through a fixed po- 
tential field V(r) remains the same if the mass of the particle 
is changed. In the corresponding quantum mechanical treatment, 
there is a direct dependence on the mass which vanishes only 

when the particle mass becomes arbitrarily large. It has been 

35 

demonstrated by Lervig et al. that electrons and positrons 
definitely exhibit many of the non classical feattrres. 
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’7.C 

Actually, it has heen estahlisbed by DeWames et al, 
that in the phenomena of channeling, the particle mass and the 
strength of the interaction potential simultaneously play an 
important role in detemining the extent to which the classical 
theory may be used. They have shown explicitly that when, either 
the particle mass or the strength of the interaction potential 
becomes small, one shoxild use quantimi mechanical treatment. This 
criteria indicates clearly as to why neutrons and protons, ins- 
pite of their equal mass, are to be treated differently. Because 
of the weak strength of interaction for neutixm case, it has to 
be treated quantum mechanically where as the protons, which inter- 
act via strong Coulomb (or screened Coulomb) field may be treated 
classically. Similarly, the electrons and positrons ara expected 
to exhibit pronounced quantum effects in the emmission experi- 
ments'^ because of their small mass. Even in the case of 
protons, where the classical treatment should be applied and the 
wave mechanical treatment has been shown"^*^ to yield classical 

results in the limit of proton channeling conditions; there are 

41 42 

some experimental observations ’ which can be explained only 

on the basis of wave mechanical theory. This point has been 

emphasized by Chadderton"^^ by considering the proton channeling 

41 

experiments of Gibbson et al. . The presence of anomalously 
high energy loss and anomalously low energy loss components in 
the transmitted energy spectra is a definite indication of the 
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Lniportaiice of interference effects due to wave nature of protons. 

Inother evidence for existence of proton waves, is observation of 
42 

star patterns . Apext from these reasons, a purely quantum me- 
chanical description of the atomic processes leading to directional 
effects in crystals is very much desirable from the standpoint of 
solid state theory. 

The wave interference in electron channeling has been treated 

by many authors'^^”'^'^ and the problem has been essentially that of 

electron diffraction. On the other hand, for the case of neutrons, 

the effects of crystal periodicity have been neglected till the 

pioneering work of DeWames et al.^^ in this direction^ published 

in 1966. They used quantum mechanics to explore the possibility 

of anomalous transmission of particles in the crystals. This has 

50 

been further investigated in detail by the author , 

In the present work, a general quantum mechanical formalism 
has been developed in Chapter II which has been applied to the 
channeling of electrons, positrons and neutrons in Chapter III. 

The problem of emission of charged particles (such as electrx>ns 
and positrons) has been treated in Chapter IV and the correspond- 
ing blocking effects investigated. The energy lost to conduc- 
tion electrons by the charged particles during the phenomena of 
channeling and blocking has been calculated using a dielectric 
formxilation in Chapter V, The conclusions have been briefly 
summarized in the last Chapter. 



CHAPTER II 


QUANTUM THEORY OE CHAJ'TNELING 

2.1 General Eoim-ul at ion; 

The Uasic concept in the theory of directional effects in 
propagation of particles through crystals is to realize the cru- 
cial role of the crystal symmetry and its effects on the parti- 
cle wave function. It has been customary in the theory of ano- 
malous transmission for particles to draw analogies from the 

corresponding effects for x-rays which is known experimentally 

SI S2 

as Boirmcari: effect • This effect is known theoretically to be 

a consequence of the crystal periodicity which leads to the for- 
mation of standing waves in the crystal when the Bmgg condition 
is satisfied. Those waves which have their antinodes at atomic 
sites are attenuated at an enhanced rate, while those having 
nodes at the atomic sites are negligibly attenuated and hence 
anomalously transmitted. The similar physical picture regarding 
the crystal periodicity is also true for the anomalous transmis- 
sion of particles in the two beam theory. Therefore every attempt 
to understand the effect theoretically, must necessarily assume 
the existence of the crystal periodicity and any deviations from 
it must be treated in perturbation approximation. This feature 

Q 

is pre^nt in the lindhard's classical theoiy (discussed in 
Sec. 1.2) in which as a first approximation, regular chains of 
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close -packed atoms jcalled strings are asstuned to exist in the 

45 

ciystal and also in the quantum formalism of DeWames et al. 
where- the ciystal is assumed to he in ground state having exact 
periodicity of lattice. 

In the present work, we' assume that the crystal is ini- 
tially in a very low lying state |n> , so that the approximate 
crystal periodicity still exists. This effectively means that 
the crystal is at very low temperature and as far as the crystal 
symmetry is concerned, the lattice vibrations do not destroy it 
to any considerable extent. Therefore, as soon as the particles 
which outside the crystal were described in terms of plane waves, 
enter into the crystal, they are influenced by the periodic inter- 
action potential of the ciystal. It is known from the band theory 
of electrons in solids that as a consequence of Bloch theorem, 
the wave function of a particle in a periodic potential should 
be a Bloch wave. Consequently, the Incident plane waves must 
go over to Bloch waves inside the crystal whose wave vectors 
should be determined by using the boundary condition at the crys- 
tal suxface. 

inother approximation that we make, is that only a few 
states which are close to the initial state of the system, are 
excited due to the incoming particles. This may not be obvious 
due to the fact that channeling takes place with fairly energetic 
particles. However, if we realize that the particles moving 
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close to the channeling conditions , do not como -very close to 
the lattice atoms, as asstmied by lindhard to illustrate the 
steering mechanism and hence they have only gentle collisions 
with the lattice atoms. Consequently the approximation is 
quite reasonable under the conditions that exist in a channel- 
ing experiment. Under this approximation, the matrix element 
of the potential between the crystal eigen states jn>and 
jm>has approximate symmetry possessed by V when jm> is closed 
to the initial state jn>and can be expanded in lattice I’ourier 
series. In particular, we will consider only one phonon inelas- 
tic processes implying thereby |m> = |n + 1 >. 

Novv we consider a perfect crystal described by the 
Schrodinger eqn. 

H^jn> = |n> (2.1) 

where is the second quantized Hamiltonian and E^ is the 
energy of the crystal wben it is in the n-th eigen state, assumed 
to be the initial state which is very low lying at some low tem- 
perature. A beam of particles each of mass m^, energy is 
incident upon it. The total Hamiltonian of the system consist- 
ing of the crystal and the particle can be written as 

H = H^ + Hp + Y (2.2) 

P P 

where H is the free particle Hamiltonian -li Y /2m and Y 
represents the interaction between the particle and the crystal. 
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Tlius the total system is described by the Schrodinger eqn. 

H 'F = Ef (2.3) 

The total wave function 'i' may be expanded as 

° m 

where r is the position of the particle, is the actual 
position of cJ -th nucleus i.e. j £ and u^ being 

the noroaeil position and the displacement of a-th atom, i«spec- 
tively. Using eqns. ( 2 . 1 )-(2,4) , we get the eqn. for ‘l>jj(r) as 

[ t v2 + (E - E^)] ♦^(r) = z (r) (2.5) 

0 m 

Eqn, (2.5) is similar to the corresponding equation obtained 
In scattering theoiy, with right hand side including a non- 
local term in n space i.e. both the elastic (m = n) and in- 
elastic (m ^ n) processes being included. Y/e will see that the 
inelastic part is responsible for both, attenuation and anomal- 
ous transmission. 

As discussed earlier, the particle wave function *^jj(3^) 
should be a Bloch wave and can be written as 

T V 

= e u^ir) (2.6) 

where the periodic function Uj^(r) has- crystal periodicity such 
that u^(r) = u^(r -i' £) , £ being a lattice vector and ^ is the 

wave vector of the particle inside the crystal which may be 
determined by using the boundary conditions at the crystal 

surface. The periodic function n^(£) can be expanded as 
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i^(r) = S uj^(n) e 

i(Kr,+l^)» r 

so that 4> (r) = z u, (n) e 

n - n 


(2.7) 


v/hare is a reciprocal lattice vector. Now, as indicated in 
the beginning of this section, assuming approximate s5njQmetiy 
for jm> close to |n> , we expand s.s 

xKw X* 

= 2 \in, m) e (2.8) 


Using equations (2.7) and (2.8) in (2.5), we get 
2 

fm -V ] ^ = ° (2-9^ 

o g,m ° 

where = 2111^ (E-l^)/ti^ = 2m^ E^ being the particle 

energy. 

Eqn. (2.9) is to be solved for the coefficients u (m). To 
do this we note that it can be written in matrix form as 

A U = - V U (2.10) 


Here U is a colimm matrix, each element being specified by the 
crystal eigen state n and the Eourier index h. Thus, a parti- 
cular element U‘ = Ujj(n) is h-th Eourier coefficient in the 
expansion of Ujj(r), the periodic part of the Bloch wave 

A A 

Similarly the elements of the matrices A and V are specified, 

A being a diagonal matrix. The partictilar matrix elements may 
be written as 


ihn 


•fi 




-4l 


^h,g °n,m 


(2.11a) 
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a^d = V, „(n,in) = fY(x) e ” ^ dr (2.11d) 

gm n-g V’ nm'— ^ ^ ^ ' 

V being the volime of the crystal. 


The ffis-trix equation (2.10) is in general, quite complic- 
cated and difficult to sol've. The high dimensionality results 
from the product of the number of phonon states excited and the 
number of Fourier coefficients to be retained in the expansions 
of and A formal solution of equation (2.10) can 

be written as 


A A 

U = Uj_ - c- V u 
.-1 


( 2 . 12 ) 


where G = (A - i£)g_^ is a diagonal matrix and is a solu- 

tion of equation (2.10) in the limit of vanishingly weak poten- 
tial, namely when the particle has not interacted v/ith the crys- 
tal and the crystal is in its initial state jn> , so that a par- 
ticular element of this column matrix can be v^ritten as 


U 


hm 




C2.13) 


The dimensionality of the matrix eqn. (2.12) can be reduced to 
a large extent by using the partitioning technique. This tech- 
nique has been successfully used in investigations of defect 
problem^^’ The idea is to partition the perturbation matrices 
such that only submatrix is significant and to keep the rest 
of the elements which are negligible, as null submatrices. Thus 
if we take only one phonon processes, the corresponding sub- 
matrix being represei)]fced as v, we can write 



y 


0 0 


(2.14) 
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Here, the matrix elements corresponding to higher phonon pro- 
cesses have been assumed to he negligible compared to one phonon 
elements contained in v. Partitioning other matrices accordingly, 
A go u 


Bi 


G U = ( ) and U = ( ,) 

0 G U’ Hi 


(2.15) 


where g, u and Uj[_ have the same dimensionality as v and corres- 
pond to the same space of n and h as v. Using these partitioned 


matrices in ( 2 . 

12 ), we get 






X. ft 

U = Ui - g V u 



( 2 . 

16a) 

and 

u* = u^ 



( 2 . 

16b) 

Now writing vu 

= s and ^ Uj_ = 

in ( 2 , 

.16a), we get 




u = u^ - g s 



( 2 . 

. 17 a) 

which yields, s 

= Mu^ with M = (I + 

A Av — 1 

^g) 

V 

( 2 . 

.17h) 


so that a particular matrix element can be written as 

„ i,h*m gm’ 

t V 2 m 

h'm _ hffl o „ 1 

& g»ni' (K^,+^) - - ie 

and using eqn. (2.13), we get for m 7 ^ n 

2 m^ u_(n) 

u, ,(m)= - p E h ' 


(2.18) 


She Bloch states ^ ^(s) excited due to propagation of particles, 

i(K^ . r 

are obtained by using (2.18) and multiplying it with e ^ 

and summing over h * , as 


2m^ 




fi 


use - >4 - le 


( 2 . 19 ) 
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Now writing eqn. (2,9) as 


2m. 




g 


h * ,mfn. 


( 2 . 20 ) 


and substituting for u^,(m) from (2.18), we get 
2 


- >4 a S Jh.g{n.n)+C(,g(r,)^Ug(n) = 0 (2.i 


o 


where 


2m. 


o 


-rrhn ,rh »m 

^h'ni “gn 


a „(n) = - Z 2 -O o ■ -- 

^ m^^n h’ "^m " 


(2.22) 


We now examine the conditions under which only a few waives 
are needed in the periodic expansion of <l>jj(r). For the ideal case, 
in which the crystal is assimied to he exactly in the ground state. 


4-9 


the conditions have been examinee, for 4^. We shall follow a 


similar procedure. If in the expansion of 4>j^(r), only first 
term u^(n) is to be significant then the solution foruj^(n) must 


be small compared to UQ(n). Thus 


Uij(2i) = ~ 




ho 


2m. 


^ * Ohh(”) 


« 1 (2.23) S 


where Uj^(n) has been written in units of u^(n). Clearly, except 


at a Bragg condition, one is dividing a Fourier coefficient of 
the interaction potential (since Clig(n) << ^ energy 


of the order of 


Eh i <S. > %-y 
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since s-s discussed below, for reciprocal lattice vectors 

perpendicular to this is about 50 eY for electrons and about 
0.025 eV for protons and neutrons. Away from the perpendicular, 
but still not satisfying a Bragg condition, for k^ >> K^, the 
second term dominates, tending towards • Thus at 

sufficiently large energy, only the planes parallel to kj^ (so 
that K^.kjj = 0) contribute appreciably to To justify the 

one wave picture, when such planes exist, it is necessary that 
the fourier coefficients of the potential be small compa3?ed 
to + h sufficiently high order reflection Y-j^ 

will in any case be small compared to Y^, limiting thereby the 
number of waves which must be taken into account. 

When the inequality (2.23) is violated for a single ref- 
lection - e.g. when a Bragg condition (hh+Snl = ’Si ^ 
filled, the first teim in the denominator of (2,23) becomes neg- 
ligible and the corresponding coefficient u^(n) becomes signifi- 
cant. This is the well known two wave picture. There will exist 
two solutions for having ^ essentially parallel to 
When such conditions are satisfied, we can write from eqn. (2.21), 
the two equations for ■UQ(n) and as 

0 

+ Y^(n)+C^^(n5u^(n) = 0 j 

- 2 . 


26 


Writing ^ + IL where n is small and its squares and higher 

powers will he neglected, we get from eqn. (2. 24) , (using the 
Bragg condition jK^+k^l = and the ^mmetry properties'^^ of 
the renormalization matrix) , 

(2.25) 

[2ilY„k^+V^(n)+Coo(n)]Uj,(n)+l¥jj(n)+0^o('')>h(”) == ° 

where n being the unit noimal to the entrance 

surface and is unit vector along k^. For nonvanishing 
values of n.^(n) and u^(n), the determinant of the equations 
(2.25) should vanish. This gives, on simplification, two 
values of n as 

= (ky2Y^Ep)5:(V^(n)+C^j^(n))-(V^(n)+C^^(n)); (2.26) 

Thus the two waves have wave vectors (k^j+H^) and (kjj+D_). The 
decay (attenuation) of these two solutions with penetration 
distance is deteimined by the imaginary part of 

Im TT^ = (k^/2Y^Ep)[+ Im Cj^^(n)-Im CQ^^(n)] (2.2?) 

We notice that for Im C^q(^) ')l T™ ^oo^^^’ solution corres- 
ponding to n will propagate almost without attenuation while 
other is attenuated at twice the rate of one wave solution. These 
solutions are, 

^ n(£) = 

s: C®^P 

^n^l) a. (1+e ”) 


( 2 . 28 ) 
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Here (r) which corresponds to ^ » has nodes at the atomic 
sites where the scattering out of the beam occurs, so that the 
attenuation is minimized, while the other wave <f>“ (r) has anti- 
notes at these sites and the attenuation is enhanced because of 
scattering. This picture is analogous to that applied in the 

cp 

theory of anomalous transmission of x-rays . It should be 
noted that even when several orders of reflection from the same 

u,- 

plane contribute appreciably to , these two solutions and 

4'” dominate and one still gets anomalous transmission. In many 

interesting cases, this two beam theory is sufficient to give 

the essential features of the phenomena. A detailed account of 

the many beam solution in electron microscopy has been given 
55 

elsewhere . 


Erom the discussion following eqn, (2,22) upto this point, 
it is clear that calculations of penetration depths and the con- 
ditions for anomalous transmission need a detailed knowledge of 
the renormalization matrix. However, one must get into numerical 
computations regarding inversion, diagonalization and multipli- 
cation of matrices to obtain further results. This expression 
(2,22) may be seen to be similar to that obtained in the corres- 
ponding exact approach to scattering theory using T-matrix 


where the problem ends up with' numerical inversion and diago- 
nalization of complicated T-matrix. Consequently, we shall now 

■■■■• ■O..... ... .. . , . . .... . . 
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2.2 First Bom-Approxi mation ; 

Let us consider the matrix M = (I + v. if the 

interaction potential is weak so that the elements of the 
matrix v are small and thu particle energy is sufficient to 
make sure that the elements of the diagonal matrix g are small 

A A 

such that the product vg may he taken to be small compared with 

A. 

the unit matrix I, then we can make a binomial expansion, to get 

A 

M = V-VgT+VgVgV- ... . (2.29) 

This expression admits of successive approximations. The use of 
the first term v yields the first Bom approximation res^^lt. 
Similarly, the second term corresponds to second Born approxi- 
mation and so on. 


Thus the first Bom approximation in (2.22) yields 
2m 


-^hn yh’m 

Chg(n) = - S S ■ ■ I? • 2 

S * m^ h' “^m 


2m 


0 1 -i(Kb+^).r+i(K^+4).r 


= - o wr /dr /dr'e 


^ 2 ( Tr I L 




^ . x 2 2 ^ 2 . 50 ) 

m?^ h» -k^-ie 


where in writing the last step, we have used the expressions for 
hn 

^h’m eqn. (2.11) and l^-(r-r’) has been suitably added 

and subtracted in the exponential. 

A 

If one retains the second term in the expansion of M, one 
gets the next order correction to Cjjg(^) as 
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(1) 2m p . 

'hg ' ^^2 ^^T’’ /ilr'/dr" e 


„ -r" 


^ ^ V (2) ''mm’(£'^ Ih (£") 

m' m^n 

i( T +kj^>( r-j ' ) ^ i( ,,+kj^) . ( r • -r '•) 

Yi' ' p p ’ ) (vT S ' ' ’ ' p o' ' ) 

^ h- (K^,+)^,)^-ic^ -le h" (K^„+k„)2-i4,-ie 


X E ® 


and in this fashion, the successive approximations may he evaluated. 
While the analytical calcula.tion becomes increasingly difficult 
for higher correction terms, the essential features may be illus- 
trated by first Bom approximation. 

As we have seen, the attenuation of the two waves in the 
two beam theory, is determined by Im n ^ which in turn needs the 
detailed calculation of imaginary part of the renormalization 
matrix elements. Irom eqn. (2.30) we write 


2m 


Im (n) = - 2 ^ /dr/dr' e 




•£+i(Kg+i^) . r* 


-fi 


. i(K. ,+]%).( r-r’) _ ^ ^ 

^ (2.31) 


Noting that the imaginary part of is anall over many unit 

cells, we put ^ in the expression (2.31). Further, since 
the summation over h’ in eqn. (2.31) extends over whole of the 
reciprocal space, we can replace it by a sum over k = Kj^i+1^ 
to write 

1 ^ ... ,2 ,.2- 


r E e 
k 


“^ml 

ik.(r-r*) ^ ^ 


( 2 . 32 ) 


30 


Since the wave vector may have any value depending upon 
the incident particle energy, the variable k is continuous 
like the usual wave vector in the reciprocal space, and one 
can replace the summation over k by an integration to obtain, 


2moTi 


Im C, (n)=--' A' /dk/dr/dr’e 




n 


2 ,2 


?, -k;.) (2-33) 

n'fn 


n ' 


.50 


which is exactly the same result as has been obtained by 
starting- with the first Bom approximation solution to 

- , I 0 1 


47c 1 r-r’ i 


%(£)=- )♦«(£’) 


in eqn. (2.5) and then using the appropriate expansions of 
and 

iTow let us choose a general interaction potential, which 
interms of the atomic potentials, may be written as V(r) = 

E representing the interaction potential 

between the incident pairticle and the atom (or ion) at position 
o. The actual position of the a -th atom R may be written 
as R(j = S + Ba » where u^ is the displacement from the equi- 
librium position of the lattice site a , due to thermal vibrau- 
tions. Thus from (2.33)j we get 


Im Cjjg(n) 


2m^Ti: 


/(% S V(k-V^h)® 


A j i z, 

Y*h^(2Ti)^ g,g’ 


To. (i-4-5g)e ' 5(k2.k2 , ) <n 1 -”l “ ' > 


X <n»|e n -g o jj3> 


(2.34) 
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where V 0 (K) = / Va(r) e^-*~ dr 

iK.u cf 

The matrix elements <n | e |n’> for n’f^n, as they are 

involved in the expression (2.34), correspond to exchange of 
phonons i.e. inelastic processes. The procedure to evaluate 
these, is to expand u^^ as. 


u .= (* )V 2 


E 

hs 


§.0 - 2. + 

C.V2 ->3 ® + aj e 3 


(2.35) 


where and f . are energy and wave vector of 3 -th phonon, res- 

tl 

pectively, and §_ ( s = 1,2,3) are their polarization vectors.The 
lattice is supposed to he monoatomic, M denoting the mass of 
each atom. There are N atoms in the ciystal. Here at and a- are 
the well known phonon creation and annihilation operators for 
3 -th phonon. 


The phonon state |n> is written as (n^,n 2 ,... >, n^,n 2 ,-.. 
representing occupation numbers corresponding to different noi®al 
modes so that the matrix element to be evaluated may be written as, 
iK. u 

-jnS = /-v* 


<n’ 


n> = <ni , ni, . . . j e 




n2,n^> 


U i(Q- a.+Qt at) 

n <n!le 3a 3 3a 3 '[n.> 
3=1 ^ 


where 


*^30 ^ 21N^’. ^ ^ ( 5 -is^® ^ 

J} s 


1/2 


if • .0 


( 2 . 36 ) 

(2.37) 


For an inelastic process involving multiphonon exchange, it has 
been shown by Kothari and Singwi'^' that the probability of any 
two phonons to have saro wave vector (i.e. exchange of any two 
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phonons corresponding to the same mode of •vibration) is negli- 
gibly small compared to the process in which all the phonons 
being exchanged, ha've different wave vectors. Therefore, for 
an 1 phonon exchange process, eqn. (2,36) inay be written as 


<n' e 


n>= n <n. + 1e jj 

3=1 ^ " 

a.+Q . at ) 




u 

n <n.|e "30 3 ^30 3' [„ > (2.38) 

3 = 1+1 ^ 3 

However, at any finite temperature T, since the exact n. corres- 

J 

ponding to a given initial state [n> is never known, we should 
take an average over initial state so that the required matrix 
element becomes 


iK.u- 
<n’ I e { n> = E 


N 


oo oo 

E ... S (I w(n.)) 

n^=0 np= 0 njj= 0 5=1 3 


-2 

1 , i(Q. a.+Q*^a'h, N , i(Q. a.+Q* at) 

n 3o 3 ^ 3 n <n^je 3cf 3 3<^ 3 jn^>(2. 


3 = 1 


3=1+1 


where m(n.) is the probability of n. phonons of wave vector f. 

0 <3 3 

being present in the field when the lattice is at a temperature 
T and is given by 


-L/k^T -n.^./k-T 
w (n.) = (1-e 3 ® ) e 3 3 B 


(2.40) 
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For an 1 phonon process this averaging has been shown^ to be 
equivalent to replacing the operator of the form exp (U) by 
(U /II ) exp (-j <U >j). Thus for emission of 1 phonons we get 


<n’ I exp ( i£.u ) ln> 

<7 5=1 


1! 


(2.41) 
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where the Dehye-Waller factor 2W is given by , 

2W = 2DK^ 

P 1 (2.42) 

with 2D = 2 ' - coth ( K ./ZtJS) 

D ^2 ^ 

We note that eqn. (2.41) gives the matrix element for 

exchange of 1 phonons each having different frequency. But, 

5Q 

since it has been shown that exchange of two phonons of same 
type is N times less probable compared to one phonon process, 
eqn. (2.41) gives the matrix element for a general 1-phonon 
process within the approximation of neglecting terms of order 
H compared to unity. 

We also see from eqn. (2.41) that an 1-phonon process 

1 /2 

involves the 1-th power of Q. -v (energy transfer/N 5 .) ' , 

Da D 

which is a small quantity, so tha,t the corresponding matrix 
element decreases rapidly as 1 increases. Moreover, the multi- 
phonon processes do not peimit even the approximate symmetiy, 
which is so crucial for anomalous transmission. Therefore as a 
first approximation, we take only the one phonon contributions 
in eqn. (2.41) so that we can write 

<n * j exp(-jZ.Ujj , ) I n> = -iQ*t (n^+1 ) ^/^ exp(-DK^) (2.43) 

where the subscript 3 = 1 has been dropped from f^ and§-|. 

The corresponding matrix element for absorption of 1- 


phonons is obtained similarly and one gets 

4 1 
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so that In one phonon approximation, we hawe 


<n| exp (iK.Up)|n’> = iQ^jn^"^^^ exp (-DK^) (2.45) 

After substituting such matrix elements in eqn. (2,34) one 
should average the expression over initial state. This is because 
the crystal is at some finite temperature T and due to thermal 
vibrations, n^, corresponding to a given initial state jn>, is 
never exactly known. Therefore we replace n^ and (n^+1) by the 
corresponding averages 


and 


<n^> = 

<n^+1>= 


5/V 


-1 


(e • - -1) 

-C/k^T -1 
(1-e ® ) 


(2.46) 


where k^ is the Boltzman constant. Thus finally we get from 
eqn. (2.34), 


Suioit 


-3 ,2^ V,(k-k^_K^)V,.(k-k^-Kg) 

,) 


V't5^(2TC)' 

2 .2 




. ^2 6(k-k^^J 

TT 1 vf ^ ^ r ^ 


‘2Mt‘ 


(1-e ® ) 


if. (2.-“') 5 ( 

'C/k-t 
(e ® -1) 


-if. (<?-<?») 
e ] 


X exp [i( k-k^ ) . (a- a ' )] (k-k^-K^^ ) . ( k-k^-Kg) 

X exp(-I)Jk-k^-K^)^ + (k-k^-Kg)^!) (2.47) 


since exp (+ = 1 . 


Since any of the phonon modes could be excited due to the 
propagation of particle through the crystal, we have summed 
over all possible values of f and divided by N, the total ni»- 
ber of possible phonon modes. The probability factors fbr tlie 
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-C/W -1 

phonon to have corresponding energy C are supplied by (1-e ) 

C/k^T -1 

and (e -1) , Thus eqn. (2.47) gives the expression for 

Im C^g(n) in one phonon approximation, for a general interaction 
potential between the particle and the crysta.1 lattice atoms and 
for a general dispersion for the lattice vibrations. 


In order to proceed a little further with eqn. (2.47)j we 
assume the crystal to have no isotopes present and that the nuc- 
lear spin of lattice ions may be neglected. Under this assump- 
tion, the Courier transform of the potential, 7^K) becomes inde- 
pendent of cr so that writing (K) = 7(K) in eqn. (2.47), we get 

2m n; 

-D-(k-^-K^)2+{k-k„-Kg)2J i(k-k^).{2-0') , ^j2 

® ^ ® T 


.2 .2 


2m^K 


a, O' 


-K"' 


) : 


2 ,2 ^“o^ 


X 


^S(k‘-k^+~^) if. 6(k^-kJ-^) 


(1-e ) 


, , , (o-o')] /p 

(e5/kB3^-1) 


4^, 


If there are Isotopes, however, additional terms will appear in 
eqn. (2.48). For example, an isotope sitting at the origin, whose 
interaction potential is SV in addition to that of a normal atom, 
will contribute terms involving integrations over the products 
of V and 6V in eqn. (2. 48). 


Now the summation over^ and a* in eqn. (2,48) can be written 
i(k -K.j:f).0, 2 

in the form js e j . This type of summations frequently 

a 

occur in the neutron scattering cross section calculations and 
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can 


be performed easily to yield ^ 2 6(lc-k +f+K, ) 


e h' 


where is volume of the unit cell. Using these simplifica- 
tions, eqn. (2,48) can he written as, 

2m % 

1“ Ohg<“) = - V.t2{2„)3 ^ 


X e 


-l‘US-i^-Kf,)^+(k-k„-K )^] , ,3 

® ■ S [(},(f) J 5(^^+f+&,) 

Ct f -u ii 

(2.49) 


c f ' * ~ h' 
^-GgCf) S 6(k-k^-Mj^,.)] 


h»t 

^2 6(k2_k2+2m^ 


where &i(f)=2M? (l.e'W) 


2 6(k2_k2-2mQ 5/h^) 


(2.50) 




Since energy C (f) is periodic in the reciprocal space, C(f)= C(f+Kjj) 
and the f-dependence of G^(f) and through ?( f ) , we 

have G^(f) = G-|(f + K^,) for all Kj. and ^2^1) “ ^2^- - ^h”^ 
all K^„. Therefore the terms in the square brackets in eqn. 

(2.49), which are being sunmed over f , can be written as 

I [|yi(f+K^,)5(k-^+f+4,)+ 2,&2(f-£j,..)6(£-£n-r'^Si,n)] 

= I S,['Jl(r«h')6(k-^+K^-^)+e2(f-K^,)6(k-Vf+rh')] 

= 2,. 'I 1 02(i-ih'>«(i-vr+4.)l 

= H I [Gi(f)6{k-^+f)-H}2(f)6(k-^-f)] 

where f is restricted to the first Brillouln zone. 


(2.51) 
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Now substituting (2,51) in eqn. (2.49), using b-fimctions, 
6(k-k^ + f) for k- integration and replacing sunsnation over f "by 
an integration V’/(2 tx)^ / df, we get 

r )Y(f+K ) 

T [ ■ (Mh)-(Mg) 


m 


1“ i ^ - [ 


2Mv ( 27 i)' 


X e 


-D f ( f-K, ) ^+( f-K„ ) o 2in^ ? 

(^hXf-Kg) e “ -S ■ 6(f2+2f.k^^ )] (2.52) 




Ibis general expression admits the use of dispersion for any 
realistic model for lattice vibrations and can be used with any 
interaction potential by using the corresponding Fourier trans- 
form in eqn. (2.52). 

2.3 Debje Model Calculations^ 

It is well known that the Debye model represents a reason- 
able approximation to the actual lattice vibrations in the crys- 
tals, particularly corresponding to the low frequency part of 

6 1 

the actual phonon dispersion. In this model, it is assumed that 
the lattice vibrates as if it were an elastic continuum but the 
frequencies of vibration cannot exceed a certain maximum value, 
chosen to make the total mamber of modes equal to the total num- 
ber of classical degrees of freedom. This is specially true for 
the long wave length lattice vibrations which cannot see the 
detailed lattice structure. If the modes are classified toy 
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their frequencies, the frequency distribution function g(w) is 

2 

proportional to w in elastic medium. We choose a maximum fre- 
quency “q a'fc which this distribution is cut off so that the 
total number of distinct modes equals 3H. The corie spending 
dispersion is simply expressed as 5(f) = = "Kef v.'here c 

is the velocity of sound in the medium. Thus the maximum value 
corresponds to a cut off in the frequency given by f^. The 
Brillouin zone, which limits the range of allowed values of f, 
is replaced by a sphere of same volme with radius f^. This 
sphere is called the Debye sphere. Since the Debye sphere is 
to conta,in N points at a density V’/S-jr in f-space, we must 
have 


i.e. 


N = 



v' ire 3 


(2.53) 


A more detailed review of the Debye theory has been given by 
Blackman^ 


In spite of the approximate nature of the Debye model, one 
gets quite satisfactory results in many interesting cases, when 
Debye model is used for lattice vibrations. Particularly, 
whenever the results depend upon an integration over f in reci- 
procal space, it is found to give very good results. The most 
striking example is that of neutron scattering cross-sections. 
The total cross-sections (elastic and inelastic »both) are found 
to agree with experiments quite well. An excellent review on 
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thermal izat ion of neutrons and applicability of Debye model 

57 

has been given by Kotbari and Singwi , 

For the present purpose we see thc,t our expression (2.52) 
contains integrations over the phonon wa.ve vector f. Therefore 
the use of the Debye model is reasonably good and is relatively 
simple to yield analytical results. This is so in view of the 
fact that particles moving near the principal ciyst allograph ic 
axes, do not come much close to the lattice atoms and very high 
frequency modes are not excited even when the particle energy 
is large. 

Now writing C = dief in eqn. (2.52)^, we get 


m df , VCf+L)T(f+K„) 

im 0^„(n)=- ° -2 “'SH-oT/kJ;' 

2Ivfficv ( 27 i)^ ^ ( 1_e ® 


X e 


6(f-2f.k^4. V ■ ^^ticf/kB!P_.V 

,2i 


(e^ 


1 ) 




-DfCf-IC )'^+(f-K^)n 2m cf 

L.-h - -g -‘6(f^+2f.k^- — ° )J (2.54) 


so that Im - 


m. 


o 


' 21®cv^(2it) 

p 2mnCf 

6(f^-2f.l^+ ^ 


/ fdf V^(f) c 


n 

-2Df^ 


and Im C 


X e 


" * (>=w:i)' ■'-* 

df Y(f+IC)¥(f) 

, (n)=- ■ - A / -? [ — 

2mQv^{2%r ^ ^ “ 

-Dr(f+5^,)2+fq^,^ 


(2.55) 


m. 


5(f2-2f.V 


X f . ( f-K|j)e 




1) 


2m cf 

V -« ■ 0 


(2.56) 
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How these expressions can be used to calculate of 

eqn. (2.27) for any interaction potential by using the corres- 
ponding Fourier transform in the expressions for Im (n) and 
Im (n) . In the next chapter we give the calculations of 
Im C^^(n) and Im ^^^(n) for neutrons and charged particles to 
find the attenuations of the two waves 4)^ and (}>~ and the con- 
ditions for anomalous transmission. 


CHAPTER III 


CAIiCULATION USING DIEEERSNT POTENTIAL MODELS 

3.1 Eermi P se udop ot ent ial - Neutrons^ 

TLe possibility of anomalous neutron transmission was 

realized shortly after this effect h;< been identified and 

51 52 

mderstood for x-rays ’ . However, no theoretical investiga- 

tion was made at that time and the only experiment directed 
toward seeing anomalous neutron transmission yielded a negative 
result , apparently due to crystal imperfections. There has been I 
no further attempt to experimentally demonstrate anomalous neutron 
transmission, The main difficulty in such experiments is that neu- | 

trons have already high penetrating power because of extremely { 

short range nuclear interaction between neutron and atomic nucleus | 
and it becomes extremely difficxilt to notice directional effects | 
in neutron transmission. However, theoretically the problem does 
carry some interest and recently, DeWames et al.'^^ have investi- 
gated the theoretical aspects of anomalous neutron penetration in 
perfect crystals. Their theory is similar to that presented in 
this work but they assumed the crystal to be initially in the 
groimd state at T = 0 and used the Einstein model for lattice 
vibrations. In what follows, w shall treat this problem using 
the formalian of Chapter II. 
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For neutrons, the two particle interaction potential is 


T?ery local, the range being of the order of nuclear dimensions. 
For wavelengths larger than the nuclear scattering length a’, 
one can replace the actual interaction hy a potential whose 
width is intermediate between the scattering length and the 
wavelength. Y/e shall choose a form which passes directly over 
to the Fermi pseudo-potential^*^ as the width vanishes: 


v(r) = 


”o 


2 2 
a e-P ^ 


(3.1) 


where a is the scattering length of a bound atom in the crystal 

and is related to the scattering length of the individual atom a’ , 

by the relation a = a' As -* 0, v(r) of eqn.(3.1) 

57 65 

goes to Fermi pseudo-potential ’ 

Thus the total interaction potential between the neutron 


and the crystal can be written as 


Y?(r) = 


2^2 2 

O 71 CT 


so that Va(r-R^) = - a^ e ^ 

o % ^ 


(3.2) 


(3.3) 


Here a cr has dimension of length and becomes scattering length 
(of a bound atom at the lattice site a ) in the case of Fermi 
pseudo-potential (p“^ -* 0). From eqn. (3.3) we calculate Ya(K) as 








a_ e 




A\ 
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In the approximation of zero nuclear spin of the ciystal atoms 
and no isotopes present, is independent of a and we get , 


Y(K) = - 


2TrB^ 


m 


a 


-£^ 4^2 . 


(3.5) 


Using this in eqn. (2.56) we get 


Im C^o(n) - 2MaQ'cT/^ ^ f 1 ® 


X 6(f -2f.k^+ 


2m cf 


2m_cf 


-D''(K.-f)^+f^] 2 

X e ^ ^ 6(f^+2f.k^- — . )] 


where 


D' = D + 1/4p' 


(3.6) 

(3.7) 


In the Dehye model, the maximum cut-off value of frequency 
f^ is always small compared to a reciprocal lattice vector K^. In 
eqn. (3.7) » we note that the integration over f will have this cut- 
off f^ as the upper limit. Consequently, if we neglect f compared 
to the reciprocal lattice vector in the Dehye-Waller factor, 
this will not cause any significant error (because the exponent 
has second power of reciprocal lattice vector). This approxima- 
tion is usually made in the calculations of neutron cross sections 
without introducing any significant error. Using this assumption 
in eqn. (3.6) we ^t, 

^3_2 df ^ (f+K. ).f 2 


57 


dicf/^T 


(3*8) 
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Now to do the angiilar integration, we choose as z-axis, 
so that = K^f(cos 0 cos a + sin 9 sin a sin 4 ) where the 

phonon wai^e wector f has been written as (1,9,4) in polar coordi- 
nates and a is the angle between and (z-axis). The coor- 
dinate system has been chosen such that is written as (Kj^,a,0), 

Thus eqn. (3.8) may be written as, 


, r^o +1 - f^+K, fxcos a . 


o ' • c -n 

f+2m c/h f^-K,fxcosa . f-2m c/h 

2 % 

where we hawe put cos 0 = x and used / sin 4 d4 = 0. Now since, for 

0 

high particle energies, is always small compared to (2k^2mQc/h), 
we can use the 6-function in angular integration to get, 

f 


,3.2, 


Tm n fy,\- Tth'^a^'e 
^ ^ho^ ^ “ 2MmQCV^ k^ 


'■ (1 + 


cos a 


2k, 


n 


) / ° f2 coth ( ) df 

0 


ffloC cos a 


2hk^ 




(3.10) 


The calcxilation of Im ^^^(n) is relatively simple. Actually 

Im C„ (n) is related to the total inelastic cross section by, 

00 ^ 


1“ - 2m;v' 


n 

2m V' 
0 


(3.11) 


where is the total inelastic cross action, ffe may calcu- 

late it directly from eqn. (2.54) without making approximation 
of neglecting f compared to in Debye-Wall^r exponentials. The 
angular integration of etpi. (2.55) yields. 
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Im C^^(n) = 


/o „2^_2D'f^ , *cfx,. /, .ON 

"T*v t J ^ (^v (3.12) 


"SMiii cv k ^ 
o c n 0 


•■2kg2)- 


This expression is more accurate than the corresponding expression 
(3.10) for Im C^Q(n) because it does not use the approximation of 
neglecting f compared to Such acciiracy is needed in calculat- 

ing the attenuation of the other wave 4> ” when 4^ is anomalously 
transmitted. 

In the low temperature limit, the f- integration can also 
be done easily, using the replacement coth (•hcf/ 2 k 2 T)^ 1 + 2 e~*^^'^^B^, 
In this limit, the De by e-Waller exponent D becomes, 






and from eqn. (3.10) we get 




Ttii a e 

2lv'Im„cv ' k^ 
o c n 


®o^^h ^ 2 

'2h k^' ^o 


3 )- 


(3.13) 


(3.14) 


with + 1/4§ . Here 9^^ is the Debye temperature defined 

as "Kcf^ = kgSj) we have neglected the terms of order higher 
than (T/Oj^)^. Fixim eqn. (3.14), we can write 

I O 

-D^ Kv, cos a 3m_c L cos a 

Im 0, .(n) = Im C (n) e ^ [l+ ■ pk , j(3. 

'00 L 2k^ 215k^f^(1-H6TVej5^)-' 


Im OQ^(n) appearing in e<pi. (3.15) is that obtained from eqn. (3. 14) 
by putting = 0 so that Bn Oj^^Cn) aad Du C^^Cn) occuring in eqo, 
(3.15) contain the same degi^ee of approximation. 
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low the wave <P ^ will experience no attenuation and will 
he anomalously transmitted when Im = 0 i.e. when Im = 

Im C^^Cn). This leads to a condition for anomalous transmission as, 


cos a 




- (1+6TVe-n^) 


3/3 = fexp (d’ _ q (3^ ^5) 


The condition (3.16) must he exactly satisfied for complete ano- 
malous penetration of 4^ and when this is so, the attenuation 
of the other wave 4” is determined by 2 Im C^qC^) ™s.y he cal- 
culated using more accurate expression (3.12). Since L.H.S. of 
eqn. (3.16) is very small for high energy particles, it will he 
satisfied only for lowest possible values of (Of course = 0 
is ideal case in which it is exactly satisfied hut physically, 

= 0 does not correspond to any plane or axis in the crystal). 

As the magnitude of increases, R.H.S. of eqn. (3.16) increases 
exponentially with exponent where as L.H.S. increases linearly 
with K^. Thus only low index (principal) crystallographic planes 
and directions may produce anomalous transmission. Moreover, 

Ip — -1 

since + 1/4p , the width of the interaction potential p 

must he anall compared to the mean s( 5 uare displacement so that 

o J 

l/4p iBay he neglected compared with D^, to write Thus 

the condition (3.16) requires that the interaction potential must 
he localized in the vicinity of lattice sites. 

When the condition (3.16) is not exactly satisfied, the 
attenuation of the beam 4^ is governed by. 
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where 


•z 2 

+ 2 _ Tct^a. 

exp _- 4 |i{j 2 ^c-f'~ y'^Ep T ( '•‘^ " 3 ) ( (3.17) 


1 2 

-DiK-u K, cos a 3m^c K. cos a 

® [1+ Sv - + ■ ■, ] 

2?ik f (1+6TVeT,^) 


(3.18) 


and the condition (3.16) corresponds to 6 = 1. As discussed ahowe, 

•j 

shorter is the range of the potential (p“ small) i.e. more loca- 
lized is thi- interaction potential, core accurately the condi- 
tion (3.16) is satisfied. In the limiting case when 1/4p may 
be neglected compared to the mean square displacement ZD^, the 
attenuation, . at short wavelength is independent of the poten- 
tial width (p~ ) and is determined only by the mean square dis- 
placement . 


The attenuation of the other wave 4)“ near the anomalous ^ 

transmission condition (3.16) is determined by 2 Im ^^^(n). This 
quantity may be calculated more accurately than Im (n) from 
eqn. (3.12). In the low temperature limit, replacing coth (-ficf/2kgT) j 
in ( 3 . 12 ), by and integrating over f, we get, I 


Im 0^„(n)=- 


r 'ilc 

*- Sk-M)'2 

10 ^ 

.2D ’ f 

1 ® « /m 


2d’ f 2 

(e 1 ° - 1 ) 


( 1+2‘e-®Il/®)+ , J,; (o^ , I erf[f (2D' ) 


■o 4D£ 


4Dj_'8D£ 


+ (2+ -fi^cVskg^T^Dj^) exp 
x(erfrf (2D’)V2+ . . tSc . l.erff ■ 


2_2/„,. 2^2 


2kg2‘(2B£) 


where erf (z) =(2/f%) f exp dt. 

0 


'2kgT<2DJ^) 


1/2 1^3 (3.19) 
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In all the abo've expressions, the limit of contact Fermi 
pseudo-potential is obtained by replacing by D^. This limit 
corresponds to maximum possibility for anomalous transmission 
condition (3.16) to be satisfied for low index planes or axes 
so th*at the wave 4 ^ is anomalously transmitted with minimum 
attenuation. On the other hand, the limit of a broad interaction 
potential (p large) is obtained by replacing by 1/4 p . But 
in this limit, the condition (3.16) can not be satisfied. This 
is because the L.H.S. of eqn. (3.16) being less than 10~ for 
typical values of and of interest, the R.H. S. cannot be 
any where close to that if 1/4 p is so large compared to that 

t 2 

£r V^-p . This is more clearly seen from Table I where (1-6) 
becomes about 10” for 1/4p = 100 Thus the condition for 

anomalous penetration requires the interaction potential to be 
short range and the atteniiation under the channeling conditions 
is determined by the mean square displaceiaent . The expression 
for e and the condition (3.16) for anomalous transmission indi- 
cate that maximum penetration occurs for low index ciystallogra- 
phic planes and directions. 

3.2 _^reened Coulomb Potential-Charged Particles: 

When the charged particles (electrons, positrons, pro- 
tons etc.) propagate through the ciystal, three types of transi- 
tions in the ciystal beccaae important: excitation of nonlocalized 
electrons, excitation of electron cloud siirrounding the lattice 
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Tatle I: Variation of (1-6) with width of the Potential 

p and temperatiire when corresponds to ^lOOl 
planes in Cu(fcc) 


r'’ 

(in units of 

jt i 


(1 - 6) in 

units of 10" 

-2 

l)^V2) 

T = 5°K 

T = 15°K 

T = 25°K 

T = 35°K 

0.0 

0.8081 

0.8145 

0.8271 

0.8461 

0. 1 

0.8101 

0.8165 

0.8292 

0.8482 

1.0 

1.0091 

1.0170 

1.0328 

1.0564 

2.0 

1.6097 

1.6223 

1.6474 

1 , 6850 

5.0 

2.6026 

2.6228 

2.6631 

2,7236 

4.0 

3.9759 

4.0065 

4.0676 

4.1593 

5.0 

5.7151 

5.7566 

5.8437 

5.9741 

6.0 

7.7957 

7.8525 

7.9698 

8.1456 

7.0 

10. 1955 

10.2693 

10.4207 

10.6474 

8.0 

12.8850 

12.9793 

13.1677 

15.4494 

9.0 

15.8380 

15.9519 

16.1795 

16.5191 

10.0 

19.0201 

19.1541 

19.4214 

19.8212 
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ion and the vihmtional transitions of the lattice ion. The first 
effect has nothing to do with the particular direction of motion 
of particle and contrihutes a large, directionally independent 
hackground to the total attenuation which is treated in Chapter V 
of the present work. The latter two effects, because of their 
localization, will shov/ anomalous behaviour at Bragg angles. In 
the present work, we illustrate the effect of localization by 
considering purely vibrational transitions, choosing screened 
Coulomb potential as the interaction between particle and crystal, 
in this section and Bom-May er potential in the next section. 


The ions of the crystal are sitting at various lattice sites. 
Their interaction with any external charged particle is expected 
to be a Coulomb potential. But because of the collective oscilla- 
tions of the conduction electron sea, this potential is screened. 
Thus the interaction between the lattice ions and the incident 
charged particle is no more a long range -Coulomb interaction 
but is a screened Coulomb potential, the strength of the screen- 
ing being determined by the density of conduction electrons. Thus 
the interaction potential between the incident charged particle 
and the crystal ions may be written as, 


2 

V(r) = z ' 

’ ^ 2 l£-lol 


so that ® 


2 e 


-A|r-R,l 


Ir-a I 

I 0 * 


= V(E) = 


( 5 . 20 ) 

( 3 . 21 ) 


and 




( 3 . 22 ) 
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Using (3.22) in eqn. (2.56), we get, 


2m^(z^Z2 ^ d£ c (f+£^).f 

T ~p^f/kjBT 


Im C^^(n) = - 


Iffli’cv. 




X 




6(f -2f.k + 


2in^cf 
O s 


n 


e 

[A^+(f+Kj^)^J (A^+f^) 


T — 

e ' 

, 2 , jj2i 


-DtCf-K^) +^] 1 

X 6(f^+2f.k - ) 7^ ' ' ■?'«?'? t(3.23) 

^ [A2+(f.K )2j(A2+f2^ j 


2in^cf 
o 


The integration inTyolved here is difficult compared to that in- 
volved in eqn. (3.7) for the neutron case. Again we neglect f 
compared to the reciprocal lattice vector in the Debyc-Waller 
exponent and also in the denominators so that the integrations 
are simplified to some extent. Again choosing k^ as z-axis and 
system axes rotated such that f = (f, 9, 4)) and = (Z^, a, 0). 
Then the angular integration may he done in a similar way as in 


the neutron case, to get finally, 


^^2 --DKh^ 




+ ( 1 + 


27 cm^(z^Z 2 e ) e 
MScVc k^ (A^+Z^^) 

Z^ cos a ^f^ f^ CQth('fecf/ 2 k 3 T) 


^ m^cZ^ cos a f^ 

■ kh k ‘ 2 ^ 

n A 


2k 


n 


■)f ° ■ -o- , 

0 (A^+f^) 


df] (3.24) 


In the low temperature limit, coth (-Rcf/kkgT) 1+26“^^^-^^®^ and 

^ 2 


f- integration may he done, to get 2 

2tim^ ( z ^ Zge ^fe "®l^h ^ c Z^ co s a 


Im (n) = -( • ^ ) 

M!5cv^ k^(A^+Zjj^) 


' 


ln(1+ ^2) 


Zj^ cos a < 5 « ^ 

+ (1 + (A , f) 1 


(3.25) 

I I i f , KANPUR^ 1 
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With T) = fQ+2fQ(T/0jj)(1-e"®DA) - A [tan"’’ 

oo . . 

where, si(x) = -/ - dt and ci(x) = -f ^4-k dt are the well 

X g? X 

known tahiilated fimctions° . In obtaining the last two terms of 
eqn. (3.26), we have extended the upper limit of f- integration 
to °o because the contribution from f greater than f^ is small, 
for the integrals in'volved there, at low temperatures. 


Nov/ from eqn. (3.25) we can write, 


-Vh" 

(n) = Im C (n) ■ -o' 2 ‘ 

°° (1+K^VA ) 


_ cos a 

i 4- 

' 2 


+(lllgOKjj cos cl/2* kjj P ^'"(A,T)) lil(1 + ° y (3.27) 

A 

so that the condition for complete anomalous transmission is. 


Kjj cos a 


2 2 2 

sc. .2 )!=[(■'+ . 2 )® -1](3.2B) 


The attenuation for the wave , which penetrates anomalously 
at the condition (3.28), is governed by, 

+ 2 - Tcm (z.Zpe^)^F^^(A,T) 

I J. ^ I ^ n . - O 1 C O ( JL r SO \ 5 _ri / ’I r\r\ > 


where 


♦hP 1 - (l-er)S-rl (3.29) 

C/ Xi p 

2 2 
-DtKv, k, cos a cos a f„ 

2 [1+ kk ■ + SO- ln(1+ -2)] (3.30) 

1+K^VA^ 2h k pf^(A,T) -* 


The effect of the screening parameter A on the attenuation 

2 

can be examined using eqn, (5. 25)* Thus for A large compared to 
so that A^» and >l> f^^, oi^ gets, 
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Im C. (n)-Ini C^^(n) = -Im C 


"Vh 


( 1 + ■ 


COS a 


y ( 3 . 31 ) 


and e®® = (1 + Kji cos exp (-D^ 

P 


(3.32) 


which is independent of ^ . Thus when the screening is weiy strong, 
implying thereby that the potential is highly localized, the condi- 
tion for anomalous transmission is independent of the screening 

A 

parameter. On the other hand if A is small compared to D-,~ and 


can be taken to be of the order of f , then, 

o ’ 


im Cjj^(n) - im 0„^(n) = -In C„^(r,) |1 - _ ^ ■ g 




1+K^ /A" 


so that 0 


SC e 


_ K, cos a m^c K, cos a In 2 

J1+ ■ V )+ ° ■ *’• • • J - (3.33) i 

2 

^ a In 2 

2 2 '(1+ 21c' ■ )+ °' " f 1 (3.3+) 

1+K^ /A - 2^1 21tt: f (1- * tan-’’ -ii) 

n o^ Tq ^ ^ 


showing that weak screening cannot lead to an anomalous effect 

SC 

because for small values of ^ of eqn.(3.34) cannot be close 

ir 

to unity and hence the condition for anomalous tienanission cannot 


be satisfied. 


3.3 Bom-Mayer Potential -Charged Particles: 

Another model potential for the interaction between the 
incident charged particle and the crystal Ions i*ich has been fre- 
quently u^d^®~'^^ as for erample in calculation of energy loss of 
charged particles during their propagation through channels, is 
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the Bora-Mayer potential. 


■VCr) = A S e 


r r-R 
c a 


(5.35) 


where represents the range of interaction. Here, both the 
constants A and p depend upon the charge numbers of the inci- 

dent particle and lattice atoms. Bor this interaction we ha^ve, 

-p i r-R I 
V (r-^) = Ae ^ 

so that Y^(K) = Y(K) = 8mA (3.36) 


Substituting this in eqn. (2.56), we get, 
SA^pZn.. df .. 


Im C, ^(m) 


0 

Who-v' 


^ ^ li!3c+(M, )2]^(?+f^)^ 


h^*— f J2. 

^ ■ r/v m 6(f^-2f.k + 

(I-s'^AbT) 

-B[(i-K.)^+f^l 

e 


2m^cf (f-4).f 


6(r+2f.k^- 


2m^cf 1 

-I” ) } 


) r (3.57) 


Again we neglect f compared to reciprocal lattice vector in the 
Debye-Waller exponentials and in the denominators, as in the case 
of screened Coulomb potential. Performing the angular integra- 
tion by choosing Ic, a.s z-axis, we get, 

41 ^ -r^rr O r\ 


Ba C^^(n) 


a,A2p2 m^e “/o' 


Kj. cos a fQ^ coth ("ficf/SkBT) 

+ (■*+ pv ) / A df 2 “>s2 ‘ ■ i 

0 (pQ )^ 


(3:38) 


In the low tempeiature limit, coth (1licf/2k^!r)j^1+e 


so that 


55 


f- integration may be performed to yield 

2 2 ^ ? 
e ^ " m^c K, cos a f 

Im Cj^ (n) = - P [; -2- 

“'=’'0 ^n(Pc«h) 2B + tj) 


+ (1 + 


cos a 


2 k 


n 


) C * ( Po >^)3 


with ( 6 

n ^ r ( 




(5.59) 


^ ,. 2 , 




+ ci( 


Pc®°v , , Pc "“ 

sinl , 


■" B " 


, --- P tic 

) sin(‘^-^) - si(-~ rr^) cos(,°m -) 


p tic 




p ^ 1 ic _ 

3 ^ 

p ^ tic _ _ p ^-! 5 c _ 




^ k T k k Ir 

B B 




Thus, Bm C (n) =- « (p^,T) 

Mhc^c ^n Pc 

and from eqn.(3.59) we can write 

2 

^-jjinv 

Im Cj,„(ii) = im 0„„{n) 


e - I > lSh ‘ 


[ 1 - 


Kj^cosa 


( 1 + K , 2 / p /)2 '■■' 


-^)]) ( 5 - 


40) 
(5.41) 


ffioC Kjj cos a f^ 
^ ' BM 




2h k^ (p .T) p/(p/+f/) 


- 


(5.42) 


The attenuation of the penetrating component 4 ) ^ is glT^en by, 

i ekp [- 

^ Y B a P 

c n p c 

2 

■DM ^>1 COS a 

where 6®^ = ^ 22 !^'*'*'' 2“'k ' ^ 

P (1+Kj^Vpc ) ^ 


m c 


COB a f ^ ^ 

211 hi ■" 


(5.44) 
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f- integration may be performed to yield 

2 2 ”^1^11^ ? 
8%A p m e m e cos a f 

Im C^Q(n) = - -9—0 o. 

p^(p^ H- fj) 


+ (1 + 


cos a 


2k 


n 


) if* (Pc’’^)] 


with J™ ^ 9 1)_ 


(3.59) 


2 /„ 2 , 


VPc+^''-^^n^/Po.^+^ /Pc) 


Po^O. -Pr'^O^ P -He 6 -He 

k^l) °°®^k^T' 


+ ciCj^j ) sUiC“-^j) - sl( 






"B-* 


pdic 

k^r 


Prv'^c PiRc BJ?5c 1 

k„T)'^°®^' k^T)'*‘®^^'k.nT)®^^^“k33? (3.40) 

JD Jd Jd 13 


Thus, Im 0 (n) =- ■ (p.,T) 


and from eqn.(3.39) we can write, 

2 


p-2lKh 

Ha Cj,„(n) = to 0„„(n) ^,2 2k. 


Kj^cosa 


n 


(3.41) 


ffloC cos a f^ 

a k„ sf*' (p^,T) p^5(p^24.f/) ^ 


c 


(3.42) 


The attenuation of the penetrating component 4) ^ is gi-ven by, 


i exp[- 




2 (1-e ®“)S.rj 


where S 


n 


BM 


“"0 ^n®p“o 
2 

-Dt K-^ rr 

^ 1 h ELu cos a 

2 ’/ if 2 'k' ) 


(3-43) 


Hi 


jC Kjj cos a f 


o ~i 


2® >S, *0 (Pc>*>Pc <Po *V> 


(3.44) 
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The condition for channeling which requires Im ^^^(n) = Im C^^Cn) 


or S = 1, is given hy 

ir 


la^c f 
o 0 


cos a _ 

■'’’"-Kp ‘(p/+£ 




wrO”,vl = hl+ -;2 O^e ■" " -f](3.4 


P 


Now we note from eqn, (3.44), that for higher particle ener- 
gies, the second and third term: may he neglected and variation of 

m'l 

the channeling parameter (1-6„ j with the range of the interaction 

Jr 

_-l 

is essentially given hy, 

d-e/”) = 1-a-®A (3.46) 

The corresponding expression (using similar approximation) for 
the case of screened Coulomh potential is obtained from eqn.(3.30) 
and is given hy, 

— i) K ^ 

( 1 _epSC) = 1 _ e ^ ^ /(1+ k// ^2) (3.47) 

The variation of (1-6 ®^) with p and that of (1-6 with A"^ 

ir ir 

at low temperatures has been given in Tables II and III respectively. 


3»4 Discussion; 

SC 

If we closely examine the expressions for ( ) and 

It 

(1-6 and the corresponding channeling conditions (3.28) and j 

P 

(3.45), we notice that the range of the interaction (p ~ and A ) 
must be small for channeling conditions to be satisfied. Appearance ^ 

p t 

of in exponential as well as multiplicative factor, at R.H.S. 
of eqn. (3.28) and (3.45) requires that the magnitude of should i 
be small i.e. the pairticles should be incident along principal 
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Table II: Variation of (1-6^*^) with width of the Screened 

P -1 

Coulomb potential A and temperature vsAien 
corresponds to {100} planes in Cu (fee). 



(1 - 


units of 10“ 

.2 

in units of 

T = 5°K 

T = 15°K 

T = 25°E 

T = 35°K 

0.1 

0.8242 

0.8307 

0.8436 

0.8629 

0,2 

0.8725 

0.8793 

0.8929 

0.9134 

0.3 

0.9528 

0.9602 

0.9751 

0.9975 

0.4 

1.065 

1.073 

1.090 

1.115 

0.5 

1.209 

1.218 

1.237 

1.266 

0.6 

1.384 

1.395 

1.417 

1.449 

0.7 

1.591 

1.603 

1.628 

1.665 

0.8 

1.828 

1,842 

1.870 

1.913 

0.9 

2.095 

2.111 

2.144 

2. 192 

1.0 

2.392 

2,411 

2.448 

2,503 

2.0 

6.855 

6.905 

7.006 

7.157 

3.0 

13.449 

13.541 

13.725 

13.999 

4.0 

21.255 

21.387 . 

21.650 

22.042 

5.0 

29.436 

29.601 • 

29.927 

30.410 

6.0 

37.388 

37.572 

37.939 

38.481 

7.0 

44.746 

44.941 

45.327 

45.896 

8.0 

51.344 

51.541 

51.930 

52.503 

9.0 

57.143 

57.336 

57.717 

58. 276 

10.0 

62. 182 

62.367 

62.652 

63.267 
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Table III: Variation of with the width of Bom 

Mayer Potential corresponding to ^0(5 

planes in Cu(fcc) 


units of 
Di) 

(1 - e 

T = 5°K 

Bl'h . 

p ) in units 

T = 15°K 

of 10"^ 

T = 25°K 

T = 35°K 

0.1 

0.8403 

0.8469 

0.8600 

0.8797 

0.2 

0.9368 

0.9441 

0.9588 

0.9807 

0.3 

1.097 

1.106 

1.123 

1.149 

0.4 

1.321 

1.332 

1.352 

1.383 

0.5 

1.608 

1.621 

1.646 

1.683 

0.6 

1.957 

1.972 

2.003 

2.048 

0.7 

2.367 

2.385 

2.422 

2.477 

0.8 

2.837 

2,869 

2.903 

2.968 

0.9 

3.365 

3.391 

3.443 

3.520 

1.0 

3.951 

5.981 

4.042 

4.132 

2.0 

12.532 

12.622 

12.799 

1 3 . 066 

3.0 

24.480 

24.635 

24.945 

25.407 

4.0 

37.487 

37.693 

38. 101 

38.707 

5.0 

49.802 

50.032 

50.488 

51.159 

6.0 

60.478 

60.708 

61.163 

61,831 

7.0 

69.221 

69.436 

69,859 

70,478 

8.0 

76,133 

76.324 

76.699 

77.247 

9.0 

81.483 

81.649 

81.973 

82, 443 

10.0 

85.581 

85.721 

85.995 

86.592 
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crystallograpliic planes or axes. These conclusions are similar 
to those derived in Section 3.1, from consideration of anomalous 
neutron transmission using Fermi pseudo-potential. A comparison 
between the results for the two potential models used for charged 
particles, shows that the screened Coulomb potential is more fa- 
vourable for channeling than the Bom-Mayer potential. This is 
.een clearly from Tables II ard III where («/=) Is waller 
than (1-6^ ) for equal values of A and showing that the 

penetration depth in the former case is greater than in later 

case. Another significant difference between the two models is ; 

1 * 

that the range A “ of screened Coulomb potential depends'only on the | 

density of conduction electrons whereas the parameters of Bom- 

Mayer interaction, A and depend upon the charge numbers of | 

both, the incident particle and the target crystal. In particular, | 

72 S 

in the Brinkman's version' of the Bom-Mayer potential, one has, ; 

A = 2.58 X 10"^ 2^72 gy 

Pc" ^eff/''-5 (3.48) 

whem a^ is the Bohr radius and z^,Z 2 are the 

charge numbers of the incident particle and the crystal atom 
respectively. Thus as z^ or Z 2 increases, decreases, result- 

ing in larger penetration depth. Similarly for small valtaes of 
z^ and Z 2 , one gets more attentiation and smaller penetration depth. 

But in the screened Goiilomb potential, there is no such dependence 
on the charge number and Zg and the penetration depth is governed 
only by the density of conduction electrons throng the scareenln# 
properties. 



CHAPTER IT 


EMISSION OP CHARGED PARTICLES PROM CRYSTALS 

4.1 Intro duction: 

It was pointed out in Sec. 3.1 tLat Lindlaard's classical 
treatment is applicable only to Heavy ion channeling. Por light 
particles (such as electrons, positrons etc.) the interference 
effects are significant, irrespective of their wavelength and 
then, one has to use the quantum theoiy developed in last two 
chapters. The esperimental situation for light -particle channel- 
ing seems to be poor. Some attempts towards such studies have 

74 

been made by Wolf et.al., who studied electron channeling in Si 

single crystals by scanning electron microscopy. Similar techni- 

73 

que has been used by Joy et al. to study the electron channeling 

in ferromagnetic crystals. Some experimental information on high 

energy electron transmission has been reported by Steeds and Taldre 

More recently, Andersen et have investigated the positron and 

electron channeling by means of Rutherford scattering. Inspite of 

78 79 

the poor experimental situation, theoretical calculations ’ 

55 

using standard many beam dynamical theory , have been made. A 

HO 

recent communication also contains such calculations for elec- 
tron channeling at high energies in gold. 

The lacE of availability of experimental data on light 
charged particle channeling would have forced one to predict 
theoretical results without comparison.' However, the diarged 
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particle emission problems which are governed by the same phy- 
sical mechanism as the channeling problems, have been interest- 
ing experimentally as well as theoretically. A detailed expe- 
rimental investigation of a-particle emission from the radioactive 

8 1 82 

nuclei embedded in the ciystals, has been carried out by Domeij ’ 

ax 

and CO workers "^5 who foimd dips in emitted intensity in low index 
directions. Gemmell and Holland have reported some measurements 
on the charged particle (proton, Deuteron and alpha particles) emis- 
sion from the crystals. These experiments have been essentially 
explained by the predictions of lindhard’s classical treatment. A 
significant work on emission of electrons and positrons from 
the crystals has been reported. by Uggerhoj''^ and coworkers Re- 
cently Walker et al.®^ have studied the effects of Bremsstrahlung 
radiations on the channeling of relativistic electrons and posi- 
trons. Although, lindhard’s classical theory has been used to 
interpret the patterns observed by Uggerhoj and Andersen ^ , 

in view of the discussion given in Sec« 1 * 3 ^ it is desirable to 
(3^0 ve lop a Q_uantuiii--]ii0 chan i cal theoiy talcing Bragg diffractions and 
the interference effects into account, DeWames et al, ’ have 
given a wave mechanical theory for this problem, but their treat- 
ment, in addition to assuming the crystal to be an Einstein crys- 
tal at zero temperature, starts with an approximate form (one 

86 

particle approximation) of the many particle Schrodinger equation 
and uses the periodic potential field of the classical model.Kiis 
one particle Schrodinger equation conserves the particle energy 
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as long as a real potential is used in the equation, and as such 
cannot he expected to include the inelastic processes and hence 
should not gi-ve rise to any attenuation. However, the attenua- 
tion of the emitted particle is incorporated by adding an imagin- 
ary part to the potential whose magnitude is chosen phcnomenolo— 
gically ♦ A detailed calculation of the attenuation “taking the 

inelastic processes into account, has been given by Humphreys 
87 

and Hirsch , Hall and Hirsch and recently by the author 
himself®^. 

In this chapter, we start with the Schrodinger equation 
for the whole system (crystal and particle) and then use the re- 
ciprocity relation to relate the problems of emission and trans- 
mission. The attenuation automatically emerges as a manifestation 
of the inelastic processes, in terms of the imaginary part of 
renormalization matrix as calculated in previous chapter. Thus 
the work presented here, provides a first principle calculation 

of the attenuation as against the earlier phenomenological treat- 

36 

ment of attenuation given by DeV/ames et al. 

4. 2 formulation : 

The process of particle emission from a crystal occurs in 
three steps: First, the particle escapes from the atom, then it 
propagates through the potential field of the crystal, and 
finally, it leaves the crystal at the surface with some probabi- 
lity distribution over the outgoing directions. Leaving aside 
the actual mechanism of emission, we shai^ be concerned only with 
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the propagation of the emitted particle through the crystal and 
its subsequent escape from the crystal. We shall also neglect 
the interactions (if any) of the emitted particle with the emitter 
itself and thus we are ignoring the processes that might occur in 
the immediate vicinity of the emitter. Under these assumptions, 
we are faced with the problem of propagation of particles through 
the periodic potential field of the crystal. The crystal is again 
assumed to be in a low lying state |n> at some low temperature, 
as in Chapter II. Therefore our system is again described by 
the Schrodinger equation (2.3). The total wave function may 
again be expanded as in eqn. (2. 4), and the equation satisfied by 
the Bloch states representing the particle wave function, is 




where = E - is the energy with which particle has been emitted 
from the crystal. 


Now to find out the intensity at a point r outside the crys- 
tal we make use of the reciprocity relation^ ’ ’ , which states 
that the intensity at a point r outside the crystal due to emission 
at a point r^ inside the crystal, is equeil to the intensity at the 
emitter site when the emitter has been placed at the observation 
peilit r. Thus 

which holds as long as one neglects the reflections at the crystal 
surface. If the R.H.S. of eqn. (4.1), representing the 
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renormalization is neglected (because its magnitude is small, as 


we will see later) the relation (4.2) follows immediately from 

^ i 

eqn. (4.1). In the general case, this relation may be derived 
by using time reversal invariance on the total 7 ?ave function 
and using the eqns. (2.1) to (2.5). The reversibility rule has 
been experimentally demonstrated^*^’^^ to hold in the emission 
and transmission problems. Therefore, the problem has now been 
reduced to that of finding the intensity at the emitter* site when 
the particles are coming from a far away point r. Ihe cLassica.1 
treatment given by lindhard, also uses this reciprocity reLation f 
to relate the two problems of emission and transmission. The prob- 
lem of particle penetration has been treated in detail in Chapters ; 
II and III. Here, we start by using the first Bom approxinatLon | 

a 

result for 4>jjj(r) ; namely, | 


2m. 


ik^ r-r' 

m 




C.r') ♦„(£') ("^•U 


on R.H.S. of eqn. (4.1) and use the usual expansions (2.6) to ( 2. 8>J 

f 

for 4>^(r) and 1^o get, i 




2 2 


where 


2m. 


-i(K. +k„).rfi(K +k ) .r' 
= - °2 ^ " 


ikjj, Ir-r' 




n*^n 


4ti| r-r' 


and Y: 


h-^ 


\_g(ii,Ja). Eqn. 


(4.4) may be written simply, as 


( 4 , 5 ) 
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(26+ ^ \-g = 0 


(4.6) 


where 26E = kf. - E, 


P 


M 


P 


(1i72m^) [K/+2Ki,.y = E 


P 


(4.7) 


and 




'h-g ~ -7-g 

Now we use the hounda.ry conditions at the surface, to calculate 
kj^. If -we restrict the sums in eqn.(4.6) to N terms, there will 
he N independent solutions for 6 and hence for the Bloch- coeffi- 
cients u-,,^; all belonging to the same energy E . These solutions 
must be superposed to obtain a wave function which matches pro- 
perly on to the plane wave which is incident on the crystal surface 
from the direction of the point of observation. Neglecting the 
reflected wave, one can write the boundary condition as, 


4> (r) = 2 A. e 


iki3)r iKv,.r _ (4.8) 




(d) e -*> - = ® 


for all r such that n.r = 0 (i.e. at the surface). One satisfies 


these botmdary conditions by choosing 


,92 


(3) = ^ +,k^ V ^.h 


and taking A- = u where the bar over 


(4.9) 

^ 3 ) denotes complex 


conjugation. In fact, the second term on R.H.S. in eqn.(4.9) is 


introduced in Chapter II, for the two beam case i.e. 


when single Bragg reflection is considtered so that 3 has only 
two values. 
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Tbe intensity at the emitter site can thus he written as, 
2 




I “■ ( 1 ) ^ ^ £ ( -i ") * — I 2 

= Is u e S u, e ^ 

j h “ 


S S e ' ^ ^ 

3 h 


>2Ile 2 S u^i^ ') 


e 


Uup>zu(3)e^^^-^] e- 

h 


r>D 

^)n.r/k ,n 


. . .(4.10) 

In the last step, the intensity contributed by the cross terms 
between different eigen solutions has been written separately in 


order to display their dependence on the distance n.r, of the 
emitter from the crystal surface. Their evaluation, while straight 
forward, is tedious and adds nothing interesting to the basic cha- 
racter of the solution. Consequently, in the usual many beam cal- 
57 

culations , a simple geometry for the ciystal and for the incident 
beam is assumed which assures that these cross terms oscillate 
rapidly with emitter depth and a collection of emitters spread 
over a small thickness will yield an average intensity which is 
just the sum of diagonal terms in eqn. (4.10). This point becomes 
more clear from the two beam theory presented below. 


The intensity pattern which one observes out side a crystal, 
due to a source placed inside, exhibits considerable structure and 
one should retain maximum: number of terms in eqn. (4.10). However, 
almost .all the features which appear, can be understood in 'teims. 
of" the ' behaviour,. .whiclj. one ..gets for .the case;bf ' 


single 


're fie ctl^ ■ = ’ 






■,,7777^" 
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Specializing to the case of an emitter located at a lattice 
site (when emitter is at ah interstitial site, the intensity maxi- 
mum hecomes minimum and ^ice versa), the two wave solution takes 
the form. 


n 


(t)l^ = x^ exp {-[1- ^ ] t/^ " } 

Vi+y2 


+ (1-x)^ exp{ -[1+ ^ 3 t/ } 

Vi+^ 


+ 2x(1-x) e 

^ 'V 

where x = ^[l+ 

fi+y^ 


-t/? ’• r ^ 

cos [v' l+y'^ t/ 


y = S/2\' > = »'h'=[V®® °ho]/®p’ 

®h = ’'b" / "o" *h = = I” ®hc/Ep) 


5 


(4.11) i 


i 

[! 

I 


(4.12) ! 

s 


and t is the distance of the emitting atom from the eiystal 
siirface. In deriving eqh. (4.11), we have made the assumption 
that = Im C|^q/Ep is small compared to which is justi- 
fied since Im is quite small compared to the Fourier coeffi- 
cient of the potential as we shall see In Sec. (4-4) . 

iEhus the calculation of |#jj(t)j^, the Intensity at the 
emitter sit© (and hence at the observation point) is essentially 
that of renormalization matrix ihad the Fourier coefficient of 
the potential Dn has been calculated in Chapter III. 
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Yie calculate now, Re using same approximations, namely 
one phonon approximation for inelastic processes and Debye 
model for lattice vibrations. 

4.3 Renormalization Matrix Elements; 

From eqn. (4.5) we can write for the real part of the 
renormalization matrix as, 

2m 

Re O^g(n) = - ^^^ 2 ° 3 f 41^ar/dr' exp[l(k-kjj-Kj,).r 

P (^2 ) ( 4 . 13 ) 


where P(1/x) represents the principal value. In writing (4.13), 
we have used the momentum representation of the free space 
Green function. 


e 



4tc I r-r ’ 


(27.)- 


/ <3^ 




For a general interaction potential V(r)=£ \ (£-%)» 

£ 

expression may be calculated using one phonon approximation for 
inelastic processes and assnming the crystal to have no isotopes 
and zero nuclear spin, as in the calculation of Im ^ 

Chapter II, we get. 


m 


CjjgCn) = - 

. Mv ( 27 i) 


" 3^ 


df [V(f+E, )¥(M„) 

T t (f+K,j).(f+Kg) 







V(f-^)V(f-Kg) 


P-2f.^+2m^£Ai^ - 1) 


-D[<I_£^)2+(f-E )2] 


+2f.lc„_2in„ ^ 


2)1(4. 


14 
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where we have snrnned over all the phonon 
the case of Im Cj^g(n). 


wave vectors f, as in 


How using the Debye model for i 

-Lattice Vibrations, we put 

5 = hcf in eqn. (4.14). Then, we use ^ « n ^ 

’ the screened Couloiiib 

potential as the interaction between th« ^ -j-x ^ , 

^ ■ci:ie emitted charged particle 

and the lattice ions as in Sec. 3.2, to get 




r ^^rU+4).(f+Kj 


TcMhcv 


p( 1 ) jS])-, 




(4.15) 


Since the maximum value of f in Debye model, f^,is small compared 
to the reciprocal lattice vector K^, we neglect f in the Debye - 
Waller exponent and in the denominators, to get 


Re C„„(ii) 


. df f.(f+K^) 

nMc(A^+K^2j J — [ 


P( , 1 ) _., 1. + , i-Cf-£h) 


f^+2f.fc - 2m of/B ^ (A®+?) 1 


(4.16) 


'rtrett: :p • = 1/-»^ Is the number density of atoms in the oiystal. 

. ^ ' V'^ r JJ 

Evaluati^^. ^he angular integration, choosing ^ as z^axis as in 
Chapter ^t. 


Ke Cjj^(n) = 


2 m 


2 s 2 -DKv ,2 


P'(z^2® ) e ^ fdf Tr „ lief 

Mftc (A +Kjj^) 0 (A'^+f^) L 


K^cosa m^cK.cosa 2k +f ^ ^ 

( 1 + 2 Y. ^..-f 2 k 3 T ^ ( 4 . 17 ) 


n 


1ik. 


n 


2 Vf 


Sow for high energy particles, where k^is large compared to f^, 
we can write In (1 + f/ 2 k^) a, f/ 2 k^ and then in the low tempera- 
ture limit coth (hcf/2k3T)2^ 1 + 2 e“^°^/^B^, one gets after integration, 

2m p'(z.z„e^)^ ^ Kv,cos a o o 09 

Re 0 (n)= ° ■ I^O+Ali )] 

me +K^ ) i 


n 


m^cKn^cos a f K^cosa p A-n/m .i '.2 

° J (1 -Ataa -1 -#)+ 2 (lH- _)f^[(,-e- 9 D /’’)(|-)2 


hk. 


n 


A--'- 2 k^~^-oL^— 


^ - 9 -n/T Ky^cosa 2 o 

■^e ^ ] -Ej^k^cos a ln( 1+ ^)+2[(1+ ) A +2K^k^cos a] 






Atic\ . ^A'SC' 


rAfiC' 


1 


( 4 . 18 ) 


Sow at the Bragg condition where the two wave picture holds, 
cos a/2k^ = energy particles, we 

neglect cos <x/2k^ compared to unity, so that eqn.( 4 . 18 ) can 


can 


be written as. 


^p'(z^Z 2 e ) 

Re Ck,^(n) = - p 2 " 

Mc(A 2 +K^ 2 )Ep 


o o 2 

2n 2^ 1 n p 2 f ^ 


7 


ln(1+ °^)] 


“h ^ 

'fc V ' , 


+ ^ e-»D/I]+ ^ K^2ln( l+f^/A^) 


' ' ‘«; 2 (i^-Kb^)[oiCj*r)ooB(^T)+si(^p)siii(^®|)] i ( 4 . 19 ) 




In the same approxliflatioa, etp. ( 3 . 25 ) for Im 
written as, ■ , ; 


can be 
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im Oj^o(n) — 


1 /? '? 0 "”^1 


Me Ep''/2(A24-K^2^ 




_A [ta.-^ ^ -l(-A^,^)sin(A«J).ai(A4°)cos(^,Jj^ 


'V' 


(4.20) 


The Fourier transform of <n|V(r)|n>, the initial state 
expectation "value of the interaction potential is gi^ven by, 


=(1A')/ (3r <n|Y(r)|n> e 

TfCll.) 


¥' 


2^ <n|e 


riK, .u 
-h -o 




n> 


(4.21) 


.. 53 


which has been calculated by Kothari and Singwi""^ and one gets, 

.2 


^h = 


47cz^Z2e 

(A^+E^^) 


p ’ e 


-W 


( 4 . 22 ) 


with W ^ 

D ^3 


S ^ coth (C./kgT) 


(4.23) 


8k temperature limit) 

B 3T D 


The equations (4.19), (4.20) and (4.22) provide the appro- 
priate values of Re Im and Y^, to be used in the expres- 

sion (4.11) for intensity. 

4.4 Discussion: 

The extinction distance ?jj’» defined in eqn. (4.12), is 
Lly less than 1(X) 1 for electron enei^ies above 100 ke? 
absorption length is quite large (more than 100 

times) td This is because Im is quite eanall 

compared below). Thus a collection of emitting atoms 



73 


spx“6a,d ovsir ss'veral ©xtinction disdaxices from ths dystad sur- 
face , such as might be achieved by bombarding the crystal with 
radio-active atoms at a definite energy, will give an intensity 
which is the average of eqn. (4.11) over thickness. This average 
will eliminate the cosine term leaving only the exponentials, 

f 

with the average distance to the emitter replacing t. 

For typical values of the parameters A = 4,2 x 10® 

= 2x1o"*® cm ^ and |V]^/Ep| = 10”^, we see on nuiKrical evalua- 
tion that the magnitude of Re is about 10“^ and that of 

^ho^^h about 10”^. Thus the attenuation is very small 
compared to what has been generally chosen phenomenologically^®. 
In the ideal case, let us study the angular characteristics of 
the intensity pattern for no attenuation (infinite absorption 
length). In this case, from eqn. (4.11), the average intensity 


is given by, 

== = x^ + (1-x)^ 

= 1 - A/d-^y^) (4.24) 

where y is approximately given by (from eqn. (4.12)), 

t' 

sin 26- 

The average intensity <j#{ > has been plotted a^tnst y in 
Fig. 1. To inc^-uie the Bragg reflection on the other side of 
9=0, the CTir^e' shotdd be' reflectsd about the valw of y. 


corresponding to - ej.ectrons, 6=0 corresponds to 

y = corresponding 


euTwe s 


"I; 

' i '.-it* 



.in figs. 2 and 3 . 


V' <h' sj;* 
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It sTaould be noted that the siiice r.e 0, ie ver/ small 

no * 

compared to these cur7es are essentially identical to those 

36 

of DeWames et al. v/itliout attenuation, Consecuentlv t^''- 
angular widths of intensity patterns are a^lso ecuaJ. to the 'widths 
of their curves corresponding to no attenuation case. However, 
their curves with attenuation correspond to an oVer estima^tion 
of the attenua,tioii 0.1) in view of the fact th»t 

Im %q/^y\ typical values chosen by le?/ames et al.^^, 


implying thereby very little attenuation. In fact, 


our 


calculation is determined once the parameters for 1', ' arc fixed 
and need not be chosen phenomenologically. The mass and energy 
dependence of the attenuation shoves that for Icr energy 

heavy particles, the attenuation becomes significant. For ke? 
electrons and positrons Re '^Yio^''^h 'patterns start 

getting modified due to this correctio'n. 

From eqn.(4.25) we note that the thiclates= averaged inten- 
sity of the two beam solution has the outer ’wing in the intensity 
pattern at (y = + 1), 




1 f ^^-"h 

32m“EYV2 

^ 0 








n 


(for electrons, at y = 1) 

2m„ 1/2 


1 r .1 . 

^ (2mQ 


(4.26) 


(4.27) 


(for positrons, at y = -1) 
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This shows that the condition A e^/e^ « 1 for vslidity of tho 

two beam theory depends upon tho mass of th^ incident particle 

and tho strength of the interaction potential V, 

^ h* 

It also depends upon the particle energy E throueh the 

P 

term Re which is im/ersely proportional to S^. The mass 
dependence shows that for heavier particle Sj as soon as A 0_ 
becomes comparable to 0^, the two beam theory is not useful. 
Similarly for strong interaction potentials large) the con- 
dition A 9^/0^ <<1 is violated. Actually the particle mass and 
strength of intera.ction potential, simultaneously govern the 
validity of the tv/o wave picture. Thus, inspite of the large 
mass, (comparable to that of proton) neutrons are describable by 
the two wave picture, since their interaction potential is v^eak 
compared to that for proton case (Coulomb or screened Coulomb). 
Similarly, the electrons and positrons, because of their small 
ma,ss5 may be described by the two beam theory (although marginally, 
because of strong interaction potential). But for protons, the 
two beam theory is certainly not apolica,ble. 

The energy dependence appears only when Re C^^^ becomes 
significant at low energies and xmder those circumstances, the 
condition A eg/e^ << 1 is more favourably satisfied for positions 
than for electrons. Another interesting consequence cf the pre- 
sent theory, in this situation, is an indication of the difference 
between widths of electron and positron emission patterais. From 
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eqn. (4.25) we see that since = (V^ + Re C^^)/Sp, it follows 
that the magnitude of y for 9 = 0 is slightly increased for posi- 
trons and slightly decreased for electrons (since Re is always 
negative). This shows that whenever Re becomes significant 
(at low energies), the angular width of electron emission pattern 
is small compared to that of positrons. Although this difference 
is only qualitative and does not quantitatively explain the obser- 

TO TQ 

ved difference"^ the difference is qualitatively correct and 

in the right direction. This feature is usually absent in the 
two beam theories^ . 



CHAPTER V 


ENERGY LOSS OP CHARGED PARTICLES 

5 . 1 Int ro du ct ion ; 

In the preceeding Chapters II, III and lY, we have consi- 
deored some of the features of the quantum theory of channeling. 
Hi particular, we examined the conditions for channeling; energj’" 
temperature and potential dependence of the phenomena and the 
widths of the emission patterns. One of the most important proh 
Lem that needs attention' in the theory of channeling, is the 
problem of energy loss of channeled particles. The quantities 
Like the actual penetration depth, the range of the channeled 
particles and the energy spectra of the emitted (or transmitted) 
particles, are essentially determined by the rate of the energy 
Loss of the channeled particles. 

Por neutral particles, such as neutrons, the only mecha- 
nism by which the channeled particles lose energj% is through 
phonon excitations. Por charged pairfcicles, the total contribu- 
tion to energy loss may be divided into three parts: (1) Energy 
loss to phonons (2) Energy loss due to core electrons and 
(3) Energy loss due to conduction electrons. 

The phonon contribution to the rate of energy loss has 
been considered by Lidiard^^. He calculated the rate of energy 
loss duB "fco Blastic processBs l)y troatiBg tbe crystal in 



81 


harmonic approximation. By regarding the fast particle coming into 
the crystal as a perturbation, it has been shora that during the 
propagation of particle, phonon shock waves are generated which 
carry away the energy from the particle. For particle velocities 
Y, much greater than the velocity of sound, the rate o'f energy’- 
loss - dE/dt was found to be inversely proportional to V confirm- 
ing the earlier results of Lehman and Leibfried^^. Moreover, the 
phonon contribution to energy loss is found to be strongly direc- 
tional effect, as expected. 

The energy loss due to core electrons is also expected to 
show directional effects, as was pointed out in Chapter III. The 
energy transferred to core electrons depends upon impact parame- 
ter and must be calculated separately for various core electrons. 
However, a comprehensive theory of this loss has not yet been 
v/orked out. It is expected that this loss should be small compa- 
red to that due to conduction electrons, because 'jnder cnanneling 

conditions, the particles do not come ver^’ close to the ions and 

93 

hence to the various core electrons. In feet Hrginsoy has shown 
that it is a good approximation to neglect this and that the energy 
loss to valence electrons is the most important contribution to 
energy loss of charged particles moving under channeliiig condi- 
tions. 

Thus the stopping power (dE/dx) of charged particles is 
mainly governed by the rate of energy loss due to conduction 
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electrons. In the first approximation these constitute a uniformly 
distributed degenerate electron gas and their contribution to the 
total stopping power is not expected to show any dependence on the 
direction of motion of the particle in the lattice. This approxi- 
mation is probably better justified in the case of metals, than 
in semiconductors. In the semiconductors, the walence electron 
density may be nonuniform so that the energy loss does depend 
upon the position of the particle. Actually, this nonuniformity 
of electron density may account for a part of the observed channel- 
ing effects. 

A general formulation of the response of an electron gas 
to a longitudinal disturbance is possible in terms of the complex 
dielectric constant e(k, w). Such a formulation has been given 
by lindhard^^ and later by ilozieres and Pines''^'"^ by using first 
order perturbation theory. The essential steps of the treatment 
given by Nozieres and Pines are as follows. 


Consider the unperturbed electron gas, described by the 
Hamiltonian, 






k k" 


k -k 


n) 


( 5 . 1 ) 


where n is the electron density and is the density fluctuations; 

given by, | 


Pk ■ ? ® 


-ik.x^ 


(5.2) 


In the equations (5.1) and (5.2), summation over i runs over all 
the electrons, whose positions are denoted by x^. . If the Fourier 
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components of the external test charge (incoming particle) are 

written as, e {r^ exp [-i( t + k.r)]+ c.c.} , then the perturted 
Hamiltonian is given hy, 


^1 , 2 ^ P _k \ ® + c. c. ) e (5. 3) 

iC 


where n is ^ infinitesimal quantity. How expanding the perturbed 
wawe function as, 

-iS^t/h 

'1' (t) = S ¥ e “ a (t) 
n 

and calculating the coefficients in first order perturbation theory 


a„(t) 


[!(-“+ “no)" ® "*1 


-0)+ - In 


-k ('’k'no ® 


W + w - in 

no “ 


(5.4) 


where (’^k^no matrix element of between the cr^J-stal 

eigen states |n> and 1 0> and ^ = (E^-E^)/h. Using the trans- 

lational invariance, which implies (p_t^)vjq = 0 if (Pk)nQ f 0> 
the perturbed expectation value of P jj- can be calculated from 
eqn. (5.4), to be (assuming that ^stem is invariant under reflection' 


<Pk> 


. r p 

*k2 ^ 


-icDt + nt 2 . 

*i<Pk)iiol ) ■ 


-) ( 




+0} -iri 
no 


low using ihe Fourier expansions of ilacroscopic Poisson eons* 

k.eE^ = 4mer^ e^“^ 

k.Ej^ = 4%e[T^ e-^“'^ + <pj^>] 
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One can write 




-icot' 


1_ 

G (kjO)) 


- I 


(5.6) 


Then from eqns. (5.5) and ( 5 . 6)5 we get, 
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Im 


__ 1 

e(k5a)')- 


4T?e‘ 

— 2“ /, 
hk"^ 


(pO 


n 


K'no^ 


(5.7) 


'jO 


Therefore, using the lowest order time dependent perturbation 
theory result for energy loss"'^^ 


one can write from eqns. ( 5 . 7 ) and '( 5 . 8 ) 

_2 


dW 

dt 


S-rce 
‘2 


Im 


1 


k‘ 


(5.9) 


6 (k,ca) 

for high velocities, where Born approximation is satisfactory each 
fourier component of incoming particle density behaves like a 
test charge of density fluctuation rj^ = 1 v/ith a frequency cj = k. V 
where Y is the particle velocity. Therefore the energj’ loss suffered 
by the fast electron is, 

,2 

k k" 


dE 

d^" 


= - 2 J®;:! (£•!) iP-eCk^fc.Y) 


1 


Replacing sunmation by integration E - ' ^ j dk, and perforraing 

k (271)^ 

angiilar integration we get, 

kV 

r ax , 

W 


dB 

d^" 


/ X e(k,'<r 

0 0 \ 7 y 


Thus the final expression for stopping power is gi^ven ty, 



dE _ 2e^ r°° dk , 1 


If the incident particle has charge z^e, the 
the integration hecomes 22^^6^/717^. 
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dw (5.11) 

factor mirltiplying 


Eor high particle velocity, the upper limit of w _ inte- 
gration may be replaced by °= and one gets - ^ a, -L . But v.hen 

ax 

the particle velocity cannot be chosen as large, the velocity 
dependence of the stopping power is modified depending on the 
form of dielectric function. Eor low particle velocities Y « Y^, 
the Fermi velocity, the problem of stopping po'ver is similar to 
that of the residual resistivity in a dilute metal alloy ?;here 
electrons are scattered by stationary charged impurities. The 
Mott^^ solution of this problem is based on the Bona approximation. 
According to lindhard^^ the same perturbation treatment is valid 
for the stopping pov/er problem because the important parameter 
in the collisions is the relative velocity. For small V, the re- 
lative velocity is simply the Fermi velocity V-^, since only the 
top most electrons can be scattered. 


The first treatment of the stopping power on this basis, 
for particle velocities small compared to the Fermi velocity, is 
due to Feossi and Teller^'^. They obtained the ejipression ^or stop- 


ping power of particle with charge ZjB^ as, 

. |£ = u -me ) (5.12) 

which shows the proportionality of the stopping pov;er to the particl 
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velocity. Here m is the electron mass, lindhard^^, Lindhard 
and Winther^® and Trulnikov amd Yavlimskii^^ have used the die- 
lectric response' loimulat ion (5.1) and derived an expression 
which differs only slightly from that of fermi and Teller, eqn. 


(5.12). For high density electron gas they obtain, 

2z 7th 

[ IJl ( 2 ) - 1 ] 


dx 


37th- 


( 5 . 13 ) 


100 


Transmission experiments in which the energy loss of protons 
and a-particles of 4-30 keV energy has been measured in thin 
films of various matrials, confirm that the stopping power is 
proportional to the particle velocity. 


In the case of heavy ions, nuclear energy losses are 

larger than the inelastic losses to electron gas, uptc energy in 

the MeV range. However, -under channeling conditions, the cross 

section of hard nuclear collisions is greatly reduced and it 

becomes possible to study the electronic stopping power down to 

relatively lowr energies ( 10-50 keV) vhere normally inelastic 

losses would be completely masked by nuclear stopping. There 

is increasing evidence that electronic excitation is appreciable 

101 102 

even for low particle velocities. Channeling experiments ’ 

with xenon ions, over a wide range of energies (1-1000 kel) in 

tungsten single crystals, show that the range of channeled ions 

changes from the E behaviour, characteristic of lew 

P C* 5 

velocities, where nuclear stopping power dominates, to the R^Ep 

behaviour, consistent with an electronic stopping povu'er proporti- 
onal to velocity, k discussion on tht velocity dependence of 



stopping power for various ranges of particle velocity, has been 
recently given by lindhard''^^. 


Here we present the results of a calculation made for 
stopping power using formula (5.11) and different models for 
dielectric function which will give us an insight, as to v/hat extent, 
the stopping power is sensitive to the fora of the dielectric func- 
tion. It is seen that for electron velocity upto about ten times 
the Permi velocity, stopping power is very nearly inversely propor- 
tional to V and is not very much sensitive to the form of dielectric 
function. It is also seen that for particle velocities small com- 
pared to Peimi velocity, even the use of expression (5.11) leads 
to the result -dE/dx very nearly proportional to velocity as v/e see 
from e qn . (5.13). 


5.2 Calc ulation Using Different Approx imat ions for Dielectric 
lection: 


for convenience of numerical computation, we ’write eqn.(5.11) 
in terms of dimensionless variables. Thus expressing lo in units of 


2 

k-p, the Permi momentum and m in units of ^k^/2m, we get, 

i? o o 


dE 

dx 


2 2 

- ' f 

2 *' 
27cK o 


dk 

k 


2kKo . 

f m im (—1.- ) 1 tu 

o £(k,ai) 


(5.14) 


writh = mY/hkj, = V/Vj,, being the ratio of pa,rticle velocity to the 
Permi velocity and Im [l/6(k,a))j is to bo 'written in the aimcnsion— 
less units mentioned above. If G(k, = £^(k, )+ iG 2 (k,w), 

e,{k,.) ana being «al, im 
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and since e 2 (k,w) vanishes for w > 2k + k^, the integral in 
eqn. (5.14) may he written as, 

2 


dB 

dz 


dk .2k+k 

L J -rr J cdim ( 


1 


2%K, 


+ / 
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£(k,aj) 


) du 


dk r 


2k K. 


2(V^) 
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J 

0 


“1“ (eCk.ojf 1 (5.15) 


Now the numerical calculation of the stopping povv'er requires the 
knowledge of dielectric function. Vfe first give different appro- 
ximations in which dielectric function has been calculated. 


(A) Random Phase Approximation (EPA) : 


The RPA was developed in the course of studies of screen- 
ing and collective behaviour of an electron liquid. In the early 
work”'*^^, primaiy emphasis was placed on the explicit introduction 
of collective coordinates (which describe the plasmcns) for this 
purpose. The approximation acquired its name on the basis of the 
physical argument that under suita.ble circumstances, a sum of 
exponential terms with randomly varying phases could be neglected 
compared to N. Thus the approximation, 


k-q 


S 

i 






N 6 


k,q 


was frequently made when appeared as a multiplicative factor 

in the basic equations of the theory. This approximation amounts 
to neglecting the exchange diagrams in the Goldstone diagram 
analysis. Thus the RPA neglects all the exchange correlations 
and represents a good approximation for high density electron 
gas (r^ < 1). 
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The expression for dielectric fimction in EPl was first 
calculated by Lindhard^"^, using a self consistent field metbod, 
This is given by, 


with , 


and 


e(k,w) = 1 + Q^(k,w) 


(5.16) 


jU+2k+k‘^ 
to -2k+k‘^ 


1 ^to-k‘^^2-| T.. ito-2k-k‘^i 
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[1-(^)U In I' 


(kjio) 


Tck' 
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71 . 


FT 


‘^2k-k‘ 


for to< 2k-k‘^ 


(5.17) 


2 2 

y [ 1-(^2k ) ] for |2k-k^l<(o < 2k+k^ (5.13) 


4k 


R/ 


= 0 


for (0 > 2k+k 


Qq (k,(o) and (k,to) representing the real and imaginary parts 


of QQ(k,a)). Here k^,^ is the Thomas-Fermi wave number in units of 
k-p and is given 

2 


, 3w 

- = = 
FT 2 2 

^ F '' F 


4 

itaokj, 


n 

2 


v/here = (47i;ne^/m) is plasma frequency. The Fermi 'wave number 

is given by''®^ k^, = (97t/4)^^^ defined as /J> = 

0 2 

r^ is generally expressed in Bohr units as r^ = T^a.^, /me 

and for metallic density region r^ lies between 1.8 and 5.6. Thu 
in EPA we have , 

(i5:,w) 

( e'(k,oj) ) ° - d'+Q S(icl)?+(Q/(k.«.))^ 


(5.19) 
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(B) Hubb^d Approximation: 

1 07 

Hubbard improved upon the results of RPA by taking into 
account the short range correlations due to exchange effects. He 
used the diagrammatic technique to include these effects by sijmming 
the exchange diagrams in an approximate way. In this approximation, 
the expression for dielectric function is obtained to be, 


with, 


e(k,oj) 

f(k) = 


= 1 + 


1-f(k) QQ(k,a)) 


2 2 
1+k"^+K 

s 


( 5 . 20 ) 

( 5 . 21 ) 


7/herc Kg is the screening parameter in units of kj,. Ibis screen- 
ing parameter did not appear in the original fomula but it v;as 
added later by Hubbard himself, in private communication to several 
authors”^^^. Although, usually the screening parameter Kg io re- 
placed by the Thomas — Fermi wave number kj^^, but it has been em- 
phasized that to reproduce the compressibility limit pro- 

perly, K should be chosen as given by Nozieres-Pines interpola- 
s 

tion formula, 


K 2 ^ _ 0^-158 k|,p)/(4 + 0.158 k|^) (5.22) 

2 

For high density metals, this formula yields vslues of Kg close 
to but somewhat less than k^^,. Thus in this approximation we get, 


Im ( 


Qq (lc,“) 


») ^ " [l4-(1-f(k))Q/(lc,<^)]^"+[(1-f(k))QQ^(k,m)] 


:5.23) 
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( C ) Kle i nman Approx imat : 

Kleinman has used self consistent method and 

later the. diagrammatic technique^ \ to calculate the dielectric 

functions talcing both exchange and Coulomb correlations into 

account. It has been shown that Hubbards formula (5.20) has 

correct form only in the static limit and that even then 

Hubbard's f(k) is completely incorrect for large k. Actually 

the static formula obtained from Kle inman 's results which has 

been confirmed by variational calculations of Langreth"’*^^, is 

found to give correct results for screening density around a 

fixed charged impurity which agree with those obta-ined by Singwi 
113 

et al. making self constent cplcula,tions. 


According to Kle inman 's formulation. 


' O V P rir 


x'^(k,-m) 

+ ) (5.24) 


^(ksW) e*(k,-w) 2 ^ e*^^(k,-m) 

'^x(k,a))+a2X*(k,-a3)+(a.|^-a2^)x(k,co) x*(k,-m) 

1+(a,j-a2)x*(k, -u) 


with = 1+ 


(5. 


where a. = 4 [l- 4 — ~] ; Up = 4[1- ^2“ ] 


(5.26) 


and X (k, +to) = x,](+w) + iX 2 (±“) where k dependence of the real 
cluantities X..|(k, +w) and x 2(^9 t)een omitted for simpli- 

city. These functions are given by, 

f oj+k^ 2 a)+2k+k^ k'^+u 1 

” u |_-^1 . i -i 


=-3 

T 2k^ ' 
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X o(") 


'-(“aV] 


if w< 2 k-k ; 0 otherv^ise 


X o(-0j) = 


Q 

! -f / ^ "*-k \ p 

‘ 2 k ^ 3 ^ -“ 2 k<(L<k + 2 k; 0 othe 


(5-28) 


rv/ise 


Here c.g^(k,w) given Hy eqn. (5.25) has significance of being 
dielectric function for an electron. This differs from the die- 
lectric function for a test charge G(k,u), given by eon. (5.24), 
since the latter has no exchange interaction with the responding 
electron gas. How using eqns. (5.24) and (5.25), we get, 


Im (l/e(lc,m)) - _(1/2) { X2(“)“X2(“‘*^) + (“-^-ci:2^ X2(“) 

-X-](‘u)X2(-^)] > +(a^-a2)^[x2(w)|x(-<^)|^-X2(-“)lx(‘^)l^I| ;/i{ 1+a 
f Xi(“) + X^(-w)]+(a^^-a2^)[x^(‘^)X^(-^)+X2(‘^)X2(-<^)] + (a^[X2(“) 

-X2(-“)^ +(a-l^-a2^nx-](-^ ) X2(‘^)-X-)(“)X2(-“)] (5.29 ) 


and Ini( 1 /eg^(k,w))= - X2(<^)-a2X2(^ )+a2(a2-a^ ) X^(‘J) X2(-w)T(a^+a2) 

2 2 2 2 — 

(a^-a2) X2(“)|x(-w)| +(2a^ -a^a2-a2 )x^ (-to ) X2( V‘ { 

’‘Ix^ (w)+Xi(-w)]+(a^^-a2^)[x^(<»j)x^(-w)+X2(‘^)X2(-“)j}^ [X2(‘^) 

-X2(~w)]+(a-i^-a2^) [ Xj(-w)X2(“)-X^(u)X2(-“)]^^i (5.3C ) 

(1) lan greth Approximation: 

Recently”^*^^, langreth has used a variational method to 
calcuaate the dielectric function. It is shown that Kleinman’s theory 
does not use the Hermitian property of interaction potentip.l. 
Consequently his results are correct only in the sta.tic case (t^i— C). 
for finite frequencies, langreth gets, 
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e(k,^^) = 1 + 




1+f(k,a)) Q^(k,co) 


(5.31) 


with f(k,<o) = 


(5.32) 


1 

where c. = -«- -■^- ; 


0 = -L 

2 2 Ttk2« 2 

s 


(5.33) 


and (k, w ), which has keen written as Xj(w) for simplicity, 


is given "by. 


XT,(k,w) 


2 fEl^'Vk^ 




(5.34) 


where f is I’ermi-function. langreth did not calculate Xjj(k/i)) 


■because he did not use it in any numerical caJ.culation. This 


may be calculated in a similar way as lindhard’s RPn calculation 


for dielectric function, The detailed steps to calculate the 


real and imaginary parts, X ^^ (‘^) and been given in 


the Appendix. The result is, 


Xt, -!(+“) 


- m [“+2“ - T ^ i 


2 1- 8k 


4 4k 


2 

Ok ) 1 


Inla3-+2k+k2|- ± [l_("^|-)Tliil“+2k+k^l for k < 2 (5.35) 


O 

i l 


kp r+^o-k 


'll Zv2 2 cjo+2k+k‘^ ^ 


4ii" 4k 


2 ^ 4k 


llnl l>for k > 2 (5.36)| 

^ w+fk+k'^ J I 


X (m) = cok^/l67tk 


for ^ < 2k-k 


w-k'" .2 


-p P w-k 

“I6xk 2k 


)^] for I 2k-k £ ^*5 5, 2k+k 


(5.37) 


= 0 


for w > 2k+k 


Xj o ( ■“ ^) 0 


(5.38) 
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Now from eqn. (5.31), we car writ( 


I+qI( 1+f^)]^ / ^ ^ i+qP( ^+pR)_Qlj,Ij 2 ^|-qI^ l+p^^)+QPpI] 2"'. (5^35^ 

where f and the real and imaginary parts of f(k,a)) in eqn. 
(5.32), are given hy, 

^ X^l(-“)]+2c2X]^^(a))xj^^(-a3)} {[Xj,^(m) 

2 2 

+ Xi,i(-“)] -Xi2(“)J H Xi2("))[xi,i(u) 

*^2 (5.40) 

P - 2 [c^X]^-](w)+Xj^2^‘*^)'^°2^I2^“^^I1^““^3 ^ [Xi’i(‘^)+Xj^^ (-i^)] 

-Xl2(“)> ^°i[Xli(“)-Xl2(“)+Xlt(-‘^)] 

+2c2Xj^^(w) x^^(-w)}: / M (5.41) 

where DN = ( [x^^^ ( ‘^)+Xj^^ (-o))] ^-X^2^ ^L 1 ^ “) 
5.3 Results and Discussion; 


Various approximations mentioned above, have been used to 
cailculate the stopping power -dE/dx of electrons, using eqn. (5. 15). 
The double integration involved in the expression has been done 
numerically, using Simpson's rule. A satisfactory convergence 
was obtained with 20 as upper limit of k- integration and interval 
0.2 for both, k and ti) , integrations. The results for electron 
velocities 2?^. and more upto those corresponding to electron 
energy about 1.5 keV, have been given in Tables IV, V, VI, VII 
and VIII, for various densities of conduction electrons (in 
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Table IV: Variation of stopping power (-dE/dx) with 

particle energy E , (r =2). 

P s 


V 


E =~ Energy Velocity 
° (eV) (loScm. 


sec""^ ) 


RPA 


(-dE/dxJ to ^ unit e of 10~^_ergs_cm”'’ 
Hubbard* HubbarS. Kle’toman Kletoinan^ 

screened 


2 

50. 12 

4.19 

14.161 

13.676 

3 

112.77 

7.97 

7.969 

7.747 

4 

200.47 

8.39 

5.114 

4.988 

5 

313.24 

10.49 

3.581 

3.500 

6 

451.06 

12.59 

2.659 

2.603 

7 

613.95 

14.68 

2.060 

2.019 

8 

801.89 

16.78 

1.648 

1.616 

9 

1014.9 

18.88 

1.351 

1.326 

10 

1253.0 

20.98 

1.125 

1.105 

11 

1516.1 

23.08 

0.936 

0.920 

12 

1804.3 

25.17 

0.787 

0.773 

13 

2117.5 

27.27 

0.671 

0.659 

14 

2455.8 

29.37 

0.578 

0. 568 

15 

2819.1 

31,47 

0.504 

0.495 

16 

3207.6 

33.57 

0.443 

0.435 

17 

3621 .0 

35.66 

0.392 

0.385 

18 

4059.6 

37.76 

0.350 

0.343 


unscre- (elec- (Test langreth 

ened tron) charge) 

13.604 S.442 12.989 

7.714 5.087 7.395 

4.970 3.187 4.782 

3.488 2.196 3.365 

2,595 1.612 2.508 

2.013 1.238 1.949 

1.612 0.984 1.563 

1.323 0.802 1.284 

1.102 0.665 1.070 

0.917 0.553 0.891 

0.771 0.464 0.749 

0.657 0,396 0.638 

0.566 0.341 0.550 

0.493 0.297 0.479 

0.434 0.261 0.421 

0.384 0.231 0.373 

0.543 0.206 0.332 


14.606 

8,226 

5.269 

3.682 

2.731 

2.113 

1.689 

1.384 

1.151 

0.958 

0.805 

0.686 

0.592 

0.515 

0.453 

0.401 

0.358 


Table V : 
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Variation of stopping power (-H/dx) v^ith 

particle energy E , (r^ = 3). 

p s 


o 

1 

t?D 

CD Cp 

Velocity_ 

EPA 

(-dE/dx) in imitsw of 10~'^ ergs 
Hubbard Hubbard Kleinman Eleinmen 

cm 





screened 

unscre- 

ened 

(elect- 

ron) 

(Test 

charge) 

Lengreth 

2 

22.28 

2.79 

8.249 

8.806 

9.474 

6.895 

9.247 

7-861 

3 

50.12 

4.19 

4.789 

5.029 

5.326 

3.651 

5.178 

4.680 

4 

89. 10 

5.59 

3.115 

3.250 

3.416 

2.269 

3.324 

3.065 

5 

139.22 

6.99 

2. 198 

2.285 

2.391 

1.551 

2.330 

2.169 

6 

200.47 

8.39 

1.642 

1.702 

1.776 

1.139 

1.732 

1.623 

7 

272.87 

9.79 

1.277 

1.321 

1.376 

0.872 

1.343 

1.264 

8 

356.39 

11.19 

1.025 

1.059 

1.100 

0.691 

1.075 

1.015 

9 

45 1.06 

12.59 

0.843 

0.869 

0,902 

0.563 

0.882 

0.835 

10 

556.87 

13.99 

0.703 

0.725 

0.751 

0.466 

0.735 

0.697 

11 

556.87 

13.99 

0.585 

0.603 

0.625 

0.387 

0.612 

0.580 

12 

801.89 

16.78 

0.492 

0.507 

0.525 

0.325 

0.514 

0.488 

13 

941.11 

18. 18 

0.419 

0.432 

0.448 

0.277 

0.438 

0.415 

14 

1091.5 

19.58 

0.361 

0.372 

0.386 

0,239 

0.378 

0.358 

15 

1253.0 

20.98 

0.314 

0.324 

0.336 

0.208 

0.329 

0.312 

16 

1425.6 

22.38 

0.277 

0.285 

0.296 

0. 183 

0.289 

0.274 

17 

1609.3 

23.78 

0.245 

0.253 

0.262 

0. 162 

0.256 

0.243 

18 

1804.3 

25.17 

0.219 

0.225 

0.233 

0.145 

0.229 

C.217 
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Table ¥1: Yariation of stopping power (-dE/dx) with particle 
energy Ep,(rg = 4). 


Z-'J Energy Velocity (-dE/dx) in units of 10"“^ ergs cn“^ 
° (eV) ( 108 cm. ■■ ■ 




see" ' ) 


Hubbard 

Hubbard Kleinman 

Ele inman 




rUrii 

screened 

unscre- 

(ele c- 

(Test 

Langreth 


. . 


— . 


ened 

tron) 

charge ) 


2 

12.53 

2. 10 

5.647 

5.984 

6. 110 

4.242 

5.849 

5.500 

3 

28. 19 

3. 15 

3.356 

3.499 

3.554 

2.320 

3.390 

3.361 

4 

50. 12 

4. 19 

2.204 

. 2.284 

2,315 

1.465 

2.214 

2.216 

5 

78.31 

5.24 

1.565 

1.615 

1.635 

1.015 

1.568 

1.576 

6 

112.77 

6.29 

1 . 172 

1.208 

1.222 

0.748 

1.174 

1.182 

7 

153.49 

7.34 

0.915 

0.941 

0.951 

0,577 

0.915 

0.922 

8 

200.47 

8.39 

0.736 

0.756 

0.763 

0.459 

0.736 

0.742 

9 

253.72 

9.44 

0.606 

0.622 

0.628 

0.375 

0. 606 

0.611 

10 

313.24 

10.49 

0.506 

0.519 

0.524 

0.311 

0. 506 

0.510 

11 

379.02 

' 11.54 

0.422 

0.432 

0.436 

0.259 

0.421 

0,425 

12 

451.06 

12.58 

0.354 

0.363 

0.367 

0.218 

0.354 

0.357 

13 

529.37 

13.64 

0.302 

0.309 

0.312 

0,185 

0.302 

0.304 

14 

613.95 

14.68 

0.260 

0.267 

0.269 

0. 160 

0. 260 

0.262 

15 

704.79 

15.73 

0.227 

0.232 

0.235 

0,139 

0.227 

C.223 

16 

801.89 

16.78 

0. 199 

0.204 

0.206 

0.122 

0. 199 

0.201 

17 

905.26 

17.83 

0.177 

0.181 

0.183 

0. 108 

0. 176 

0. 178 

18 

1014.9 

18,88 

0.157 

0. 161 

0. 163 

0.097 

0.157 

0. 159 

19 

1130.8 

19.93 

0. 141 

0.145 

0. 146 

0.087 

0. 141 

0. 142 

20. 

1252.9 

20.98 

0.128 

0. 131 

0.132 

0.078 

0. 127 

0. 129 

21 

1381.4 

22.03 

0.116 

0.119 

0. 120 

0.071 

0. 116 

0.117 

22 

1516. 1 

23.08 

0.105 

0.108 

0. 109 

0.065 

0. 105 

0. 106 

23 

1657.0 

24. 13 

0.096 

0.099 

0. 100 

0.059 

0.096 

0.097 

24 

1804.2 

25.17 

0.089 

0.091 

0.092 

0.0 5't 

0.086 

0.089 
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tTable VIII: Variation of stopping power (-dS/dx) 





with 

particle 

energy 1 


6 ). 



^o" 

■ J - Energy 

Velocity 


(-dS/dx) 

in -unii: 

;s of 10 “ 

ergs cn' 

-1 



''e (eV) 

( I0°cm. 
se c- ”* ) 

REA 

Hubhard 

screened 

Hubbard 

unscre- 

Kleinman 

(elec- 

, Kleinnan 
(Test 

Lang ret 







ened 

tron) 

charge ) 



2 

5.569 

1.39 

3 . 

422 

3.546 

3-625 

2.490 

3.449 

3 . 

103 

3 

12.53 

2.10 

2 . 

091 

2.140 

2.173 

1.389 

2.045 

2 . 

033 

4 

22.28 

2.79 

1 . 

388 

1.414 

1-433 

0.887 

1.352 

1 . 

367 

5 

34.80 

3.49 

0 . 

991 

1 . 008 

1.020 

0.619 

0.965 

0 . 

980 

6 

50. 12 

4. 19 

0 . 

746 

0.757 

0 . 766 

0.453 

0.726 

0 . 

740 

7 

68.22 

4.89 

0. 

583 

0.592 

0.598 

0 . 355 

0.569 

0 . 

579 

8 

89. 10 

5.59 

0 . 

470 

0.477 

0.481 

0.283 

0.459 

0 . 

467 

9 

112.77 

6.29 

0 . 

388 

0.393 

0.397 

0.232 

0.379 

0 . 

386 

10 

139.22 

6.99 

0 . 

324 

0.329 

0.332 

0. 193 

0.317 

0 . 

323 

11 

168.45 

7.69 

0. 

270 

0.274 

0.276 

0 . 160 

0.264 

0. 

269 

12 

200.47 

8.39 

0 . 

227 

0.230 

0.232 

0.135 

0.222 

0 . 

226 

13 

235.28 

9.09 

0 . 

19 ^ i - 

0.196 

0. 198 

0.115 

0. 189 

0 . 

193 

14 

272.87 

9.79 

0 . 

167 

0.169 

0. 171 

0.099 

0. 163 

0 . 

166 

15 

313.24 

10.49 

0 . 

145 

0. 147 

0 . U 9 

0.086 

0 . K 2 

0 . 

145 

16 

356.40 

11 . 19 

0 . 

128 

0.129 

0.131 

0.076 

0. 125 

0 . 

127 

17 

402 • 34 

11.89 

0 . 

113 

0.115 

0.116 

0.067 

0.111 

0 . 

113 

18 

451.06 

12.59 

0. 

101 

0. 102 

0. 103 

0.060 

0.099 

0 . 

, 100 

19 

502.58 

13.29 

0 . 

091 

0.092 

0.093 

0.054 

0.088 

n 

•U t 

,090 

20 

556.87 

13.99 

0 . 

.082 

0.083 

0.064 

0.049 

C . C 8 C 

U i 

,081 

21 

613.95 

14.68 

0 , 

.075 

0.075 

0.076 

0.044 

0.072 

0 . 

.074 

22 

673.81 

15.38 

0 . 

.068 

0.068 

0.069 

0.040 

C .066 

0 , 

.067 

23 

736.46 

16.08 

0 . 

.062 

0.063 

0.063 

0 . C 37 

0.060 

n 

.062 

2 ^ 

801.89 

16.78 

0 

.057 

0.058 

0.058 

0 . 03 - 

0.055 

'O 

.057 



1C0 


me-tallic density range) i.e. for r = 2 to r = 6, .i eorparison 

leetween the results of different approximations has been made 

for a particular electron density (r^ = 3) in Pig. 4 v;he-re _dE/dx 

has been plotted against The gr^ph shows that results of all 

the approximations for dielectric function shov; essentially the 

same behaviour i.e. the stopping power to be very nearly inversely 

proportional to the particle velocity in the velocity range of 

2Yj, to about 10 Vp, Por higher velocities, the ■';6lccity dependence 

tends to become 1/V^, as expected, [in Pig. 4? curve 6 corresponds 

to use of eqn. (5.30) in eqn. (5.14)’]. These results establish 

that the velocity dependence of the energj" loss due to conduction 

electrons is rather independent of the model chosen for dielectric 

function. The magnitudes of stopping power also do net differ 

cuppreciably when different dielectric functions are used. This 

result is to be contrasted with the fact that static impurity 

screening is highly sensitive to the particular model used for 

1 1 ?_ 1 14 

dielectric function where only static dielectric constant 

(ci3= o) is required. 

finally, we have given few results of stopping power of keV 
xenon ions iri tungsten targets in Ta.ble IP. The velocity of the 
xenon ions is veiy small compOired to the Peimi velocity. According- 
ly , -dE/dx is found tc be very nearly proportional xe the particle 
velocity (shov;n in Pig. 5) as also expected from the formula 
(3.12) or (3.13). 
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Table IX: Variation of stopping power (-dE/dx) of 
Xe in tungsten ciysta.1. 



Energy 

(keV) 


Velocity 

(10® cm. 
sec“’' ) 


(-dS/dx) in units of 

. -4 -1 

10 ergs_ cm 

RPA Hubbard (screened) 


0.010 

0.467 

0.026 

0. 198 

0.226 

0.012 

0.675 

0.051 

0.255 

0.295 

0.014 

0.916 

0.056 

0.288 

0.551 

0.016 

1.196 

0.041 

0.525 

0.575 

0.018 

1.514 

0.047 

0.555 

C . 408 

0.020 

1.869 

0.052 

0.565 

0.421 

0.022 

2.262 

0.057 

0.425 

0.438 

0.024 

2.691 

0.062 

' j . 46 5 

0.554 





CHAPTER VI 


CONCLUSION 


Several interesting conclusions regarding the plienorr^ene- 
of channeling of particles in perfect crystals, emerge fron 
the results of the investigations presented in this work. The 
fact that the present quantum-mechanical theory confirms some 
of the results expected classically, is in itself interesting 
and is in the usual spirit of quantum mechanics. Thus the 
validity of the condition for anomalous transmission, requires 
that the magnitude of the reciproca,! lattice vector corres- 
ponding to the plane (or axis) v/ith which the iz'.cident particles 
are being diffracted under the Bragg condition | i = 1 | » 

should he smallest possible. Por energetic particle (in which 
case is large compared to K^) this raeans that the incident 
particles should enter into the crystal very nearly parallel to 
the principal (or lov/ index) crystallographic directions or 
planes. This is what v;a,s expected classically in view of the 
fact that low index planes end strings are separated by maximum 
distance between consecutive planes or strings end hence give- 
rise to open channels through which pairticles may propagate with 
minimum attenuation. This aspect is implicit in the Lindhard’s 
classical treatment. The effects of lattice vibrations are 
found to decrease the channeling effects and the temperature 
dependence appears as Debye-Waller factor, as one should have expe 
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There are se^ral significart results obtained in the pre- 
sent forrnalism which change many of the cla.ssical notions, alto- 
gether. It is usually argued that for high energy particles, 
since the deBroglie wavelength of the particle is sra.ll corpared 
to the lattice spacing, there should not be any significant inter- 
ference effects. On the other hand it is well known that the 
phenomena observed in the electron microscope must be described 
in terms of wave interference: Bragg angles and resonance widths 
dominate the intensity patterns .even when the electron v/avelength 
is very small compared to the lattice spacing of the target crys- 
tal. To this end, the present quantum-mechanicel treatment shcv/s ' 
that the ezistence of anomalous effect has nothing to do with the 
wavelength of the incident pa,rticle, but rather is determined by 
two criteria, namely that the interaction potential must be suffi- 
ciently weak compared to the particle energy and that the inter- 
action be localized in the vicinity of the lattice sites. If the ' 
range of the interaction is small compared to the m.ean square dis-^ 
placement, as is expected to be c=.se for neutrons, the a.ttenua- 
tion of both the waves and ^ ^ bear.: theory, is 

independent of the potential width and is determined only by the ■ 
mean square displacement. On the other hand, for long range 
interaction potentials, the condition for anomalous transmission 
cannot bo satisfied. Comparison bet?t?een the results of two 
potential models used for charged particles, the screened Coulomb 
and the Bom-Mayer, shows that the former is more favourable for 
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channeling. Even for equal range, the penetration depth is 
larger for screened. Ooulomh potential than for the Boin-L'ayer 
potential. It is also seen that the channeling condition has 
some energy dependent tc-ms. Although these tertrs are snail at 
high particle energies, they become significant to reduce the 
channeling effect, at low pa,rticl6 energies. Si'nlarly, in 
addition to the usual Debye-baller factor, a snail temperature 
dependence comes in the channeling-parar.iete-r (1-0) which does 
not show any significant effect at low t empe natures. 


It is interesting to note that energy dependence (1/Ep) 
of the exponents in the equations (3.17), (5.29) and (3.43), for 
hhsn the channeling parax'nter (1-0) is not exactly sere 


but a small quantity i.e. channeling condition is approximately 


satisfied, the wave 4(1 will channel the particles but there 
¥/ill be some aftenua.t ion. Iho attenuation uptc a given aistmce 
traversed will be governed by dist-nce at -which 

a given fraction of the initially channelea particles is present 
in the channel, is proportional tc E^. Ihis conclusion is in 
confirmation with the e 2 i)ertmentally esaablisnea fact uha.v. 
the thickness into the ciystal at which one half of the initially 
channeled particles have escaped from the channel is proportional , 
to the particle energy E^. 

The quantum theory of the directional effects in the 
charged particle emission in the crystals, presented in the 
Chapter IV, also leads to very importajit conclusions. The 
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mass and the potential dependence of the angles 9^. ^nd A in 

ij J:5 

equations (4.26) and (4.27) shows that whenever tho particle 
mass or the strength of potentirul hocomes small, the interfe- 
rence snd Bragg diffraction effects modify the emitted intensity 
pattern significantly (since A becomes small and the 

classical orbital picture is no more applicable.) This com- 
pletely modifies the classical notion of particle trajectory 
being independent of mass for given energy and potent ia.1. The 
classical limit is regained only v/hen the particle mass becomes 
very la,rge compared to the electron mass (sore than 2C0 to 300 
times electron mass) because in this limit A 0^ becomes compa- 
rable to ©2 and two beam theory loses its validity. Ilorecver, 
in this formula,tion, thu attenu'^tion of the emitted particle 
emerges automa,tica,lly a,s a iricnifestation of inel^^stic processes 
in terms of imaginary part of renormalization r.atrix and need 
not be chosen phenomenologically. The actual =^ttenuaticn is 
foimd to be sm.all compared to what is usua.lly chosen phencmenolc- 
gic.ally and is proportional to (m^/E^)"*/^. The t emperat^ore de- 
pendence of the a.ttcnuation appears a,s Dabyo-V/aller factor. 
Another interesting feature of the present formalism is an indi- 
cation of the difference in the wadths of electron ■:'rd positron 
emission patterns even in the tv/o be^n theory. The magnitude of 
this difference is dete mined by the magnitude of Re which 
is relatively small but has ezi energy' dependence . Thus at 

low particle energies, one can see tho-t the positron emission 
pattern has more width as compared to that of electron emission 
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pattei?!!, as has "been experi-nentally observed^®’ Although, 
the present results do not quant itavely explain the experrrienta.lly 
obsorvod difference, they are certainly in the right direction 
and give a qualitative explanation v.'hich is usually absent in 
the earlier tvro beam theories^® ’ 

The co,lculation of energy'' loss of charged particles pre- 
sented in Chapter V shows that the contribution of conduction 
electrons to the stopping power is inversely proportional to 
the pc^rticle velocity for velocities greater than the ronii ve- 
locity upto about 10 tines the I'erni velocity/. For higher pc.rti- 

0 

cle velocities, the dependence goes over to 1/7^. Fron Fig. 4-, 
it is also clear that the velccitj^ dependence of the stopping 
power is rather insensitive to the particular approxinaticn used 
for the dielectric function. The actual magnitude is also not 
very sensitive to the fonn of dielectric function. 

Thus we see that the present quantur. nechanical treatnent 
of channeling phenomena in crystcals leads to roany new interesting 
fea.tures, net expected classically. It ‘is expected that a nore 
detailed analysis using nany beer theory ray yield few rore fine 
structures, with the same basic characteristics obtained here 
in the tW'O bear theory. More realistic rxdels for lattice vib- 
rations ray also be examined using formulation of Chapter II. 

The problem of energy loss of channeled particles requires a 
little more attention. Firstly an elegant formulation for the 
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contribution of coro electrons to energy less is very nuch 
desirable because it is expected to shov directional effects, 
in addition to giving a correction to energy loss duc. to con- 
duction electrons. Secondly', as bas been pointed out recently 
that even the conduction electron contribution is also expected 
to show some directional effects and stopping power nay be 
slightly different imder chcnneling conditions than in normal 
conditions. This is significant in sejuiconductors where the 
density of conduction electrons is not uniforr.. Therefore, it 
is desirable to improve upon the dielectric femui^tion of 
energy loss due to conduction electrons and to develop a unified 
formulation for the energj'’ loss of charged particle s due to- con- 
duction and core electrons. 
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APPEiJDlX 


If we write (o , p and k in energy units (i.e. = -p^ etc.)> 


the function Xj^(k,w) may be written as 


Xj^(k, w) = 
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how using the property of Fermi- function f„ at zero temperature, 

where f = 1 for p < k^ and 0 othervilse, we can v.'rite from (A.2), 
p u 

the real and imaginary parts as, 
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* Here -x^ = (k^,^ - P^ - k^)/2pk. 
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Xl'i(“‘^) is obtained by replacing w by - in tbe e5q5ressions 
(A. 3), (A. 4) and (A. 5J while becomes zero because of 

the product of Permi functions ip( ^®sn used in angular 
integration with +2p.k) , gives zero. 

Now performing the angular integrations in eqns. (A. 3) to 
(A. 6), one gets, 
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The p- integrations may also be done easily and finally we get. 
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and X " f^Tik ' for“< 2kkj>-k^ 

2 

= 16^ J 2kVfe^<“<2BCp+k2 

j5 

If we write these equations in the dimensionless \ 
used in the text, we get eqns. (5.35) to (5.38). 
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